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ABSTRACT 


Historically, the inventory control problem has been 8 
analyzed under the assumption that the exact quantity or- 
dered was delivered in one lot, at the end of a probabilistic 
or deterministic lead time,or, alternatively, that receipts 
are the direct result of a continuous, fixed rate production 
process. In this paper, the effect of relaxation of these 
assumptions is considered. It is quite common in large in- 
ventory systems to have the total quantity ordered deliv- 
ered in increments phased over time. Contracts providing 
for delivery of short or excess quantities within tolerances 
are also frequently encountered. Production rates may be 
stochastic due to work stoppages, the diversion of resources 
to meet higher priority requirements and other external fac- ¥ 
tors. A family of models, with a variety of delivery con- 
ditions, are presented. The optimal operating policies and 


costs of these models are compared to the policies and costs 


resulting from models with conventional delivery assumptions. 
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CHAPTER 1 
INTRODUCTION 

1.1 Background 

The inventory control problem is as old as history 
itself, but the first analytic studies of significance were 
those of Ford Harris and R. H. Wilson about 1915. These 
simple lot size models determine an optimal inventory pol- 
icy in terms of a single variable, "The Economic Order 
Quantity." The assumptions required for this model are 
quite strict, and yet this simple model is useful in many 
applications. It is an excellent model for uncomplicated 
inventory systems where there is little motivation for 
assuming stochastic variation in demand. Efforts between 
world wars extended this simple deterministic model but it 
was not until the growth of the Management Sciences and 
Operations Research after World War II that attention was 
given to the stochastic nature of the problem and to the 
complex nature of many inventory systems. Recently multi- 
echelon systems have been examined and stochastic demand 
and lead time models have been analyzed. Consideration 
has been given to the special management problems of re- 
pairable item systems and other complicated systems. It 
is perhaps surprising, then, that virtually all of the 
literature deals with models for which very strict deliv- 
ery assumptions are made. 


It is almost universally assumed that either the 


exact quantity ordered arrives at a single instant at the 





end of a (stochastic or deterministic) lead time, or, 
alternatively, that the exact quantity ordered becomes 
available at a continuous fixed rate, starting at the end 
of a (stochastic or deterministic) lead time. Note that the 
first assumption is a special case of the second, with the 
delivery rate infinite. It has also been generally assumed 
that an order is placed with a single source, selected for 
offering the lowest unit price or by some other procedure 
outside of the model. Little attention has been given to 
determining an optimal operating policy under other deliv- 
ery assumptions. Existing literature which does apply to 
this problem will be referenced when appropriate. 
1.2 Methodology and Assumptions 

The conventional delivery assumption actually consists 
of several ideas. First, the quantity is assumed to be 
exactly identical to the quantity ordered. Second, the 
receipt rate is assumed to be known, constant and con- 
tinuous with the infinite rate being the most usual version. 
Third, delivery is assumed to be from one source, generally 
chosen outside of the model. Relaxation of each aspect of 
this assumption will be discussed, although this paper 
does not presume to remove all delivery restrictions or 
even to totally remove any one aspect of the restriction. 

The basic approach will be to consider models which : 
have relatively strict assumptions concerning other fea- 


tures of the system. This approach is necessary; other- 


wise the models would be so complex as to make the effects 





of the revised delivery assumptions very difficult to sep- 
arate from other unrelated effects. In general, the models 
presented are either exact deterministic models or approx- 
imate probabilistic models. The occasional approximate 
deterministic model will be clearly labeled as to the 
approximations used and the reasons for doing so. In some 
cases, the probabilistic models are so designed as to pro- 
vide qualitative information only, and some very rough 
heuristic rules are given as the only quantitative results. 
In some areas this type of result is very useful and exact 
models would have to assume knowledge of distributions and 
parameters which are not usually known. 

The basic technique will be to determine optimal 
values for the decision variables, where the optimal values 
are taken to be those which minimize the average variable 
cost. This is equivalent to minimizing the average total 
cost. The phrase, decision variable, is used here to de- 
scribe those parameters over which the inventory manager 
has control. The setting of values for these parameters 
is the way in which an inventory manager exerts his in- 


fluence over the inventory system. A specific set of values 





for these parameters, in terms of the fixed inputs, con- 
stitutes an operating policy. The average costs, referred 
to above, are time averages over a sufficient period of 
time to cover several order cycles. In this paper, except 


where specifically discussed, it will always be assumed 


that the fixed input parameters will remain constant over 





a sufficiently long period of time to allow this assump- 
tion: time averages are an accurate reflection of the 
costs to be minimized. 

Throughout this paper, decision variables will often 
be considered as if they were continuous when, in fact, 
they obviously take on only discrete integer values. The 
models could, of course, be formulated for integer solu- 
tion, using difference equations. However, continuous 
models are generally used because it is easier to deter- 
mine optimal values using the calculus. It is pointed out 
that the general problem addressed in this paper makes 
little sense unless the variables have fairly large values. 
Unless the order quantity and the demand rate are fairly 
large numbers, there is little reason to consider any 
delivery pattern but the usual one. The systems for which r 
this paper's results are most significant and applicable 
will be ones for which demand is hundreds or thousands 
per year and order quantities are of comparable size. 

Lead ‘times are not specifically considered for the 
deterministic models. While this is not an assumption in 
the conventional sense, it is more conveniently considered 
at this point. The lead time,t, is defined as the time 
lapse between the placement of an order and the receipt of 
the first increment of material. In deterministic demand 
and lead time models, knowing the point in time at which the 


first receipt is required, it is easy to determine the time 


at which an order should be placed. Equivalently, knowing 





the amount of stock on hand required at the time of first 
receipt allows computation of a reorder level in terms of 
either inventory position (IP), or net inventory (NI). 

The basic factors which allow the lead times to be assumed 
equal to zero in these models are: that IP and NI differ 
only be the on order (00) quantity which has a time aver- 
age of At, the lead time demand, and that each cycle will 
be identical to every other cycle. It is then possible to 
determine t’, such that t' = t-mT, and mT < t < (m+l1)T, 
where m is an integer. Now t' lies between zero and T. 

We can then compute areorder time or reorder point using 
t'. When either demand or lead time is stochastic, de- 
tailed consideration is required within the model. 

One assumption is not usually considered in the lit- 
erature, but this author's experience with military in- 
ventory management prompts its inclusion. The fiscal 
climate must be assumed adequate for optimal or at least 
constrained optimal operation. Whether external con- 
straints such as warehouse space, total inventory value or 
number of contracts per year are imposed or not, a minimal 
cost will be derived, subject to the conditions imposed. 

It must be assumed that the controller of the purse strings 

will provide the money necessary to cover these costs or, 

when appropriate, will close down the system. This assump- 

tion is basically equivalent to assuming that the system is | 


operated as a commercial enterprize or, at least, is using 


profit maximization or cost minimization as the basic goal, 





It is fairly obvious that inventory systems are not 
necessarily orientated primarily toward fiscal profit and 
loss. Military inventory systems are, of necessity, pri- 
marily oriented toward effective support of the operating 
forces in war or peace, with cost minimization as a sec- 
ondary goal. A much broader approach is required to 
evaluate optimal policy in a trade off between maximal 
support (supply availability) and minimal cost, since the 
former is not easily translated into dollar terms. It 
is also obvious that certain items, because of their ex- 
tremely critical nature, will be held at high inventory 
levels, regardless of the economic cost in the narrow 
sense. The climate in which the models of this paper 
operate will be assumed to be subject to constraint only 
in such limited senses as floor space or on hand value. 
A constraint on dollars spent will not be allowed, although 
such constraints exist in reality. 
1.3 Conventional Delivery Models 

For the convenience of the reader, the following re- 
sults from inventory theory are given since they will be 
referenced often. They are the "simple lot size" or "Wil- 
son Q" model and a simple deterministic finite production 
rate model. 

The Wilson economic order quantity model minimizes H 
average total cost for a system where no shortages are 
allowed and the entire quantity ordered is delivered at 


a Single point in time. Demand is deterministic. It is 


Le 





assumed that the reader is familiar with this basic model 
and the purpose of presentation is to provide easy refer- 
ence in the notation to be used throughout this paper and 
to introduce the notation in a familiar context. Figure 
1 is a graphic representation of the model. Since the on 
order quantity is constant, it need not be considered in 


the cost minimization. The only decision variable is Q. 


net Ke Total cost per cycle, where a cycle is the 

time period between successive orders or 
equivalently between successive deliveries. 

K = Variable cost per unit time. 

Q = Order quantity. 

c = Unit price. 

A = Order cost. 

r ='Demand rate (deterministic). 

z = Holding cost rate (applies to on hand only). 

b -= Buffer stock (on hand at time of initial 
receipt). 

i = As superscript implies an optimal value. 


The total cost is the sum of order cost, holding cost 
and the cost of the items. Now 


KIe =A+ rc(g + b)T +0C. 


It is obvious, by inspection, that the cost will be minimal 
when b = 0, since b < 0 would imply shortages exist which 
is prohibited. To find K it is necessary to divide xo” 


by T = Q/\A and drop the fixed and zero cost terms. The 


variable cost per unit time is 
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FIGURE 1 


WILSON "Q" MODEL 





To find the optimal value of the decision variable Q, set the 


derivative of K with respect to Q equal to zero and solve for 


©. 
ak ai”. ie 
ee aoe, Be nih 
d “eh 
ar ig 
implies 
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 \ ic * 
d*K _ 2ad 
——* =o. °° 
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implies this value of Q minimizes K. Now the optimal values 


are seen to be 


ozo = (SR, 
m= 0, 
and 
KY = k= VeaA Le % 


The simple deterministic finite production rate model 
varies from the Wilson Q model in that delivery is assumed to 
be the result of a continuous production process which allows 
immediate placement into stock. Figure 2 illustrates this 
model. Inventory level builds up at a rate of ¥ - A during 
the production portion of the cycle and then is depleted at 
the rate } during the remainder of the cycle. Two new sym- 
bols must be defined. Let ¥ = continuous production rate as 
defined earlier and P = length of time during which produc- 
tion takes place. Now 


x = oC + AY ICT E SA 4 b| ; 


LS 
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FIGURE 2 


FINITE PRODUCTION MODEL 








Again, it is obvious that b* = 0. We divide by T = Q/d and 
drop the fixed and zero cost terms to find the variable cost 


per unit time, 


x = BA, LcQ cS 


Set the derivative with respect to Q equal to zero and solve 
for Q. The second derivative with respect to Q is always 
greater than or equal to zero, so this value of Q yields a 


Minimum K. The optimal values are 


b®¥ = 0» 5 
2? = Or FY , 
and l 
y 2 
% 2 = 
Ke <= Ro K a] 
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CHAPTER 2 

INTRODUCTION TO PHASED DELIVERY MODEL STRUCTURE 
2.1 General Discussion and Definitions i 

In chapters two through five, several models are pre- 
sented. The assumption of a continuous finite or infinite 
delivery rate is relaxed, to allow incremental deliveries 
phased over a period of time. It is assumed that period- 
ically, an incremental quantity Qs will be received and after 
an interval T; another quantity Asay will be received. The 
leadtime is Tos Recall the lead time was previously defined 
to be the time between placement of an order and the receipt 
of the first delivery of units. The total quantity ordered, 
Q, will be assumed to be delivered in n lots of size qe such 


that 


after the order is placed. The delivery is phased over a 
cycle time of 
TT. 


n 
o tT. , 


where Th is the interval between the delivery of the nth (last) 
Lee qn and the placement of the next order. 

We will first consider models for which the delivery lot . 
size and time intervals between them are equal for all values 


of i. That is q, =4q for i = 1,2,...,n and Ms = t for 


L8 





y 


Sed ya ST ott. We will then crmmuntet models in which the Wore 
sizes and intervals may be aibiierarg.. | 

One common assumption is made for all of the following 
phased delivery models. When the number of deliveries is a 
(potential) variable, the order cost is of the torm A + nB 
where A is order cost occurring as a result of placing the 
order and B is the component of order cost resulting from the 
delivery of an increment or lot. Teo a wexy close approxima-— 
tion the elements of order cost over a period of time are 
either proportional to the number of —— placed ee may 
be included in A, ot are proportional to the be oe of ship- 
ments received and nay be included in B, Or are proportional 
to the number of units purchased and may be included _ the 
unit cost C of the item. For example, the administrative 
cost of letting a contract and set-up costs are included in 
A, the rim ostertiwe cost of shipping and receiving a lot 
is included in B and the carrier's transportation charge ver 
pound is included in C. 

Inventory systems often receive ordered quantities 1n 
increments and could do so more frequently it sonsamced 4 
was desirable. " Many inventory managers have recognized the 
value of placing a wontract against which calls may be made, 
which is a type of phased delivery. In fact, on g May 195 7., 
Cyrus Vance, then Deputy Secretary of the Department of De- 
fense, ina memorandum to the Secretary of the Army stated 


that deliveries from procurement should be scheduled into 


stock in proportion to anticipated issues, since, in this 





way, Operating levels may be controlled without effecting the 
total order quantity or frequency of reordering. 

A manufacturer may have some or all of the quantity of 
the item ordered on hand in his own inventory system, or he 
may have to produce the item after receiving the order. There 
is, in general, an administrative lead: time-and :potentially 
(unless the entire quantity is available in’ the producer's 
inventory) a production lead’ time until =the ‘first ‘item | ftexcept 
for a partial quantity on hand) comes off the end of an 
assembly line completed. All items may be completed as a 
batch, but more frequently, assembly line production will im- 
ply a finite discrete or continuous production rate. We will 
define a discrete delivery rate R which is dependent upon the 
production rate ¥. However, R need not be equal to ¥ and may 
be either larger or smaller than ¥. Usually, in current 
practice, shipment is made in one batch after the entire 
quantity has been produced leading to R =~, the case of a 
single delivery which is covered well in the literature. 

The assumption that the delivery rate is continuous and 
finite has also been well covered, but this assumption is 
unrealistic unless deliveries occur at very short time inter- 
vals (one a day perhaps). Charges for shipment must be in- 
significant or dependent only on the total units shipped. A 
combination warehouse factory or a warehouse adjoining a fac- 
tory meets these criteria and the finite rate model is prob- 
ably excellent for use in managing a producer's inventory of 
finished products located in a factory warehouse. However, 


this assumption rapidly loses its applicability, as distance 


20 
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enters the picture, or whenever two or more distinct activ- 
ities are involved. If the source is some distance from the 
inventory location, it will not usually be feasible to ship 
continuously and we must consider R, the discrete analogue 

- of ¥. This paper will, hereafter, reserve the notation R for 
a discrete delivery rate and use ¥ when a continuous produc- 
tion and delivery procedure is considered or as a limiting 
value for R. Basically, the reader is expected to think of 
Ras fixed size shipments arriving at fixed time intervals. 
For example, R = 200 per month implies one shipment of 200 
units is received once each month. 

In addition, phased deliveries will usually decrease 
lead time which is sitions universally accepted as a desired 
objective in "real world" systems. Obviously, it will be pos- 
sible in many cases to have material shipped as it is produced 
rather than waiting for all units to be ready for <a s 
This has the advantage of effectively reducing the lead time 
for a portion of the order quantity. 

The models which follow represent a selection of theoret- 
ical and practical systems. The first is primarily designed 
as a starting point for analysis and to show the relationship 
between phased delivery and the conventional delivery assump- 
tion. The later models are designed to represent real sys- 
tems and give practical guidance for setting operating policy 
under the conditions assumed. 

2.2 A Simple Theoretical Model 
This model is a no shortage model which assumes the man- 


ager has the option of getting some number, n, of equal sized 


aL 


shipments, q, at equal intervals, t. The effect of the man- 
ager's control over each of the parameters (number, size and 
interval) is investigated. 

The decision variables are n, q and t. Shipment will 
be made in n equal sized lots, each of q units, at intervals 
tT. In the conventional delivery models n =1, q = Q and tT 
is undefined. Note that At > q is infeasible since the demand 
rate \ exceeds the supply rate R = q/t. Therefore it will be 
required that q > At. Then Q = qn and T = qn/d\. The order 


cost per cycle is A + Bn, so the order cost per unit time is 


(At+Bn)) _ AX , BA 
qn qv gy 


Holding costs per cycle are proportional to the area under 
the net inventory curve. Reference to Figure 3 shows that 

an easy way to get this area is by subtracting the small 
parallelograms from the larger triangle with area ; QT. Each 


parallelogram has area tq and there are 


of them. Thus holding cost per cycle is 


i= 


n-1 
IC S -qt J | = =$(or-qtn(n-1)) 
1 


= = (qnT - TAt(n-1)), . 


and holding cost per unit time is 


5 (qn - At(n-1)). 





QUANTITY 











FIGURE 3 


N-INCREMENT DELIVERY MODEL 
(N=5) 





Combining these terms the average variable cost is: 


= ae Br iC@n _ ICni\t IC\t 
oe? eS ee ee 


If n 2? 1 then K is minimal when t is maximum. But \T < q, 
therefore, t = q/A when n > l. 

When n = 1 the formulation reduces to the Wilson Q model, 
q=Q=Q, with an order cost A' = A+B, The interval be- 
tween deliveries, tT, is indeterminate when n= 1. (Define 
t = 0 when n = 1 for convenience of notation.) 

However, we are primarily interested in the case when 
n # 1 and K is minimum for tT = q/\. That is tT* = q*/d. The 
problem has now been reduced to two decision variables since 
t* has been shown to be a function of q*. 

Now we can rewrite K since when Tt = q/A, 


aes (qn-ATnN+AT) = = , 


so that K becomes 


Now set the derivatives of K with respect to gq andn 


equal to zero and solve. 


— «= + === 0 
P) 2 2 
| wa 4 
so that 1 
2 (A+nB)]2 
gq [2tasne) . (2.1) 
oK _ -AA 
ie = OF 4 (2.2) 4 
nq 


The Hessian is positive definite so we have a minimum. Now 


equation (2.2) implies that either n* = ~ and q* > 0 or else 


24 





q* = © and n* > 0. But n cannot be less than 1 so from equa- 


tion (2.1) 


2 _ 2AA , 2BA 
? sac. Ic 


< 2X (asB) . 

Therefore q is finite and the optimum must be with n* = ~, 
Now n* = ~ implies l 
The optimal value of K is found by substituting n* and q* 


into the equation for kK. 


i 2.2...2 2 i 
Ke AATEC +2 BY} iby i — 2BA)? 
me A Lk pea Jae Wek Bel : 
oL2BA 2Bir 2*1I¢c 
or 1 
K = (2B 1c) 7 ; 


This result is of the same form as the Wilson Q, but has B 

in place of A. This is not unexpected, since A is spent only 
once while B is spent as a result of each delivery. 0 — 

is infinite and it seems quite reasonable that the Wilson 
formula apply to each increment rather than to the infinite 
order quantity, Q. However, such an infinite order quantity 
is obviously not applicable to any real system. Indeed, this 
model has ignored many realistic constraints. Conditions vary 
over time and organizations are reluctant to commit themselves 
over long periods. It is Phe ford important to consider the 
results when n is constant (not a decision variable) or con- 
strained. Other constraints will be considered in a later 
model. 


When n is constrained to some maximum value N, the optimal 


solution from equation (2.2) is for nq to be a maximum which 


20 








Occurs when n is equal to its maximum allowed value N. If n 


is constant (n=N) or restricted to finite values, (n<N) then 


al 
2A (A+NB) 2 
~ Moe oo 4 
n* =N ¢ 
and 1 


A! _ [2N) (A+NB) on 
ey ao [awe 


The optimal value for K is found by substituting n* and q* 


into the equation for K. 


1 1 eh 
ae a*,*IcN og peer | 2 [z-co2a (aren) 12 
N22) (A+NB) 2) (A+NB) NIC z 
or dL 
Ke = [AS ee | 2 
~ N 


This value of K* is always less than the conventional Ky 
since (A+NB)/N is’ less than A' = A + B, the order cost for 
the Wilson Q model. 


If B = O then 


Peale ’ 
N= 
5 
a — 
Qe ange: , 
and 
K,, 
Ke = 
Ne 


These results appear qualitatively reasonable since holding 
cost is reduced relative to the Wilson Q formulation by the 


incremental delivery schedule, while with B = 0,order cost is 


unchanged. It can be seen that this will increase total order 





quantity but require q be less than Q., and result in an 


optimal cost, K*, which is less than Ky the Wilson cost. 
2.3 Comparison with Finite Production Rate Model 

It is also informative to note the relationship of 
this model to the finite production rate model. Consider 
the optimal n given fixed values for q and t. Recall 
R = q/t is the discrete analogue of ¥ and that R-¥ as Tt and 
q>0O and n>». Also recall 


Ad B He 
AN, BAY = (qn=—nAT+AT) . 


K= 
nq q 


If we let B equal zero then 


= <a . 3c 4 
se = a + =H (q-At) = 0. 
n“q 


Solving for n we find that 


2Ar 
ke = 
Po [- Teste | 


since 


implies that a minimum exists at the solution point. Now 


Q = nq so that 


i ea lee Goan . [3a aS! . . 


Note that this can be written as 


Nie 
=) 


v= 0, (eta) 


This result is consistent with the finite production model 


with R replacing ¥. It is exactly the same result in the 


limit as Tt and q go to zero and n goes to infinity. 








CHAPTER 3 
CONSTRAINED PHASED DELIVERY MODELS WITH EQUAL 
DELIVERY LOT SIZE AT EQUAL INTERVALS 
3.1 General Discussion 

It is obvious from the preceding model that an uncon- 
strained incremental delivery model is applicable only 
when the manager has very extensive control over the system 
and knowledge of the present and future values of the param- 
eters involved. Although the inventory manager will very 
seldom have such complete control of the delivery process, 
he will usually have some degree of control or a set of 
choices he can exercise. 

The following phased delivery models assume varying 
degrees of control, and apply realistic constraints to the 
system. The models have increased management control in 
order of presentation. The first section actually contains 
two models. We first assume Tt and q are fixed constants 
set by the manufacturer (only n is variable) and then assume 
that the production (delivery) rate q/t is fixed as a ratio 
but that subject to q/t equals R we can vary q, T, andn. 
This is a restricted 2 variable model. The next model in 
the series fixes only the interval t between deliveries and 
leaves q and n as variables. The final model is a three 
variable model with Tt, q and n restricted only by inequal- 
ity constraints. 


These models are deterministic models designed for 


application wherever expected values can be used in lieu of 





random variables for demand, leadtime, and order quantity. 
Shortages will not be allowed in these models. Probabil- 

istic considerations and a discussion of shortages will be 
considered in sections 5.1 and 5.2. 

The models are similar to the model of section 2.2, 
but the effect of the varying degree of management control 
and the constraints is to limit the search for optimal 
operating policies to an area of feasible solutions. This 
helps to insure applicable, implementable results. 

Some discussion of desirable constraints is neces- 
sary. The most obvious constraints are that n and q must 
be positive integers. For the reasons discussed earlier, 
we will not require n and gq to be integers but the con- 
straints n = 1 and q 2 1 prove useful. A constraint on T 
of the form tT 2 4 puts an upper bound on delivery frequen- 
cy which may be necessary to prevent an unreasonable work- 
load on the receiving department. Direct constraints on q 
of the form gq s M may be necessary to prevent overloading 
of warehouses or to keep the shelf time down to prevent 
Spollage. Constraints of the form q 2 mM, may be neces- 
sary to insure carload lot size shipping economies or to 
keep the value above some minimum order size. For simpli- 
fication, recall that q 2 ar, sot 2 a implies 4T 2,4a=m 


2 


which is equivalent to q =m If we have several con- 


straints of the form q > m., we can combine them into one 


constraint q 2 m where m is the maximum of the m,;. In 


29 


tI (ha A 


this way constraints on t and q can be consolidated into a 
pair of constraints q >mandq < M. 

In a real system uncertainty always exists and there 
is always an explicit or implied planning horizon which 
may not be exceeded. In this context, a constraint on 
maximum time for use of ordered assets is imposed. This 
constraint may be applied in the form Q/A < T or ng < AT. 
Obviously the time constraint is exactly equivalent to a 


constraint on maximum order quantity nq = Q applied di- 


: r ; ; 
rectly. (i.e. Q < 7 = Wy) A constraint on maximum order 


value would be exactly equivalent in form. If the maximum 


order quantity must be less than X dollars, QC < X is 
exactly equivalent to Q < n = Wo. If several such con- 


straints were to be applied, they could all be put into 


the form Q s Ws and the single constraint, Q < W = min (W, ) 
i 


insures all are met. This constraint,Q = nq < Wiis non- 
linear in terms of the decision variables, n:>and q. Sim- 
ilarly a nonlinear constraint of the form w < O== ng could 
be imposed to keep down the workload on the purchasing 
department. 

An upper bound on n, which was applied earlier, is 
rather artificial, since the constraints on t and Q prob- 


ably represent the true motivation for constraining n, 


Note that the constraints used imply 1 <n<s-. 
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A constraint on maximum on hand quantity would be a 
constraint of the form n(q-At) < z. This constraint is 
not imposed here to prevent overcomplicating the model, but 
it is obvious that the constraints imposed define an upper 
bound on the maximum stock on hand. Since ng = Q < W, 
n 21 and) > 0 it is obvious that for n = 1 the maximum 
on hand is Q < W, immediately after receipt of the quantity 
Q. If n> 1 then t > O and the maximum on hand quantity 
will be n(q-At) = Q - nAt < W- Xt << W. The maximum on 
hand quantity occurs immediately after the nth receipt. 
When n is greater than one, the maximum on hand will prob- 
able be very much less than W. 
3.2 Delivery Size and Interval Fixed 

Suppose the source of material offers shipments of 200 
on the first of each month. Now q and Tt are given constant 
parameters and our only decision variable is n. Recall 
that for feasibility, q must be greater or equal to AT. 
See Figure 4 for a graphic representation of a cycle. 

From the prior model we have total variable cost 


ye ee ee 
= oe! ag eS (nq-nAT+AT). 


Since n is the only decision variable, the variable cost K 


can be rewritten as: 
K = —= + = (q-\At)n 


To find the optimal n, we set the derivative of K with 


respect to n equal to zero and solve for n, checking to 
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FIGURE 4 


DELIVERY SIZE AND INTERVAL FIXED 
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insure the second derivative of K with respect to n is 
positive to insure a minimum. Now, 


oK -Ad IC 


ia ie Se 2 
nq 
implies 1 
at = [eo 2M i" 
ICq (q=AT) E 
Also 
a°K _ +2A0A , g 
2 3. 
an nq 


for ali nm < ™ 

There will be no problem with imaginary roots since 
q > At and all parameters are strictly positive. If 
n' < 1,then n* = 1 and t = 0, otherwise n* = n'. 


If n* = n' 2 1,then substituting n* in the equation 


for K we get 
1 1 
Bee 


Zz 2 2 = 
= (HX Tq Wea5) BA (IC) “(q-At) “2Ad CAT 
eS a DAA a” [See | *—T* 


or 


{e 


2 
K* =[zarze US: ia * et) | ~ a + aes - 


which may be greater or less than K_,, depending upon the 
values of the input parameters. The manager has insuffi- 
cient control over the situation to guarantee a cost less 
than Ko: 


tf n" < i then a* = 2, t = © ana 


_ (A+B)A mK a 
. a + =sh > kK 


and is equal to K, if and only if q is equal to o., 2 
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Before passing on to the next point, it is worth 


noting that if q = At + e€ then lim n* =o . Also 
e>0 


_ fBARICE , BA IC (qte) 
K* = ‘a + a + =e 7 


and 


Bid Be 
Lan 3 eee = 
e+0 q ‘ 


which is of the form of KX and is minimal when 


| flax 
\ 1a * 


Now let us suppose that the manufacturing source offers 
shipments at any intervals but will ship all items pro- 
duced since the last shipment. That is q = Yt where ¥ is 
the production rate. q is a decision variable but t is 
restricted by the relationship q = ¥T, so we can substi- 
tute g/¥ for t For feasibility ¥ 2 4. Then 


FB Ped o PO hag at 
K = a + = ia 3 (ng nay + vy): 


Now we take the derivatives of K with respect to n and q 


and set them equal to zero to solve for the optimal values. 


oK _ -AA =nBA 
a as 2 2 
nq nq 





TG jg TAA Rs ues 
+ = (hn->- + 7) = 0 
implies 


oo 2) (A+nB) 
_ TAA 
nic (n->-+7) 


oK _ AA IGG 7,_A) « 
Pal caea + = tag) Ae 
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per af ee ee etn. Ley 


fa eee hw gen 


see weenie. sale pe hetero 4 


f 


implies 


n2 2Ai ( ¥ ) 
Yeas 


or 


- am ea 
I¢Cn 2 


2° an =A 
The Hessian is always positive definite, so that we are 
assured of having a minimum. Setting the two equations for 
q? in terms of n equal to each other, we get n* in terms of 


the invariant parameters, which can be substituted into 


either equation to find q*. We now have 


and 


which imply 


i 
2 
= [yee 
‘i 3 i i te 7 
It is probably not surprising to see that 
iL. 
2 


nam = [egal 


which is the same result as for the finite production rate 
model of the conventional literature (see Section 1.3). 

This means that we can decide on our total order quantity 
in the usual fashion and then ex post facto decide on our 


delivery schedule using the preceding model. 
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However, it is probably most surprising to find q*, 
the optimal delivery lot size, is a function of ¥, the pro- 
duction rate, and not X, the demand rate. It has the usual 
form for an optimal order size formula otherwise. The i 
dependency is built into the optimal value of n*. Perhavs 
the best intuitive explanation for this is that the choice 
of q* and t* is forced on us by the delivery condition and 
our demand rate must exert its influence through the value 
of n. 

In a draft Army Material Command report, Deemer and 
Hoekstra [3] consider a probabilistic version of this model. 
They simplify the model by assuming that 1t is always fixed 
at one month and that delivery quantity is not only constant, 
but exactly equal to AT, the expected monthly demand. The 
probabilistic nature of demand implies an increasing risk of 
stock out, since the expected safety stock remains constant 
throughout the contract period. This results in a high pro- 
tection level during the early part of the contract period, 
and lower protection late in the cycle. It would appear 
better to choose a delivery lot size somewhat in excess of 
At. This will result in the expected value of the safety 
level increasing as the demand variability increases. For 
a given expenditure for holding safety stock, we would be 
able to maintain a fixed protection level, or make an approp- 
riate trade off between rising risk and rising stock levels. 

For example, consider a one year contract period with 


monthly deliveries. Assume demand is normally distributed 
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with a time dependent standard deviation o = t, where t is 
in months. We decide to hold 156 unit months of safety stock. 
If we establish a constant safety level of 13 units with a 
equal to the expected monthly demand, then the average pro- 
tection level is 
1/12 (1.0 + 1.0 + 1.0 + .9994 + .9953 + .9849 + .9684 
+ .9479 + .9257 + .9032 + .8814 + . 8606) 

= .9639 
If we establish an increasing safety level by starting with 
a safety level of two and let q be equal to the expected 
monthly demand plus two units, the safety level during the 
ith month will be 2i. Since the standard deviation in the 
ith month is i, the average protection level will be constant 
at .9773. 

In this example, we can increase the expected perform- 
ance with the same average safety level by using an order 
quantity greater than the monthly demand and allowing the 
safety level to build up during the contract period. 

3.3 Delivery interval fixed 

In this model it is assumed that Tt is fixed, while n 
and q are decision variables. This type of problem would 
arise whenever shipping schedules were externally set. Ac- 
tually there may be a choice of several values for Tt. For 
example, a heavy item may be transported by ship with a twice 
a month schedule. Then our choices for t are ¥ month, 1 
month, 1.% months, etc. An item which is batch produced may 
have a 4 month batch processing time. If only one batch can 


be produced at a time, we could, under certain conditions, 
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model this with t equals 4 months, 8 months, etc. In some 
applications, management policy may specify a t which is 
based on factors exogenous to the model. Top management may 
specify monthly deliveries of raw materials based on a com- 
bination of shipping schedules, storage capacity, and the ex- 
tra workload due to more frequent deliveries. If more than 
one value of tT is possible, the solution procedure must be 
repeated for each value of Tt and the proper value of t chosen 
by ‘comparison of the costs of the resulting optimal operating 
policies. It will be seen later that only part of the compu- 
tation need be recalculated. 

Recall gq must be greater or equal to At for feasibility 
and 


ek: ee 
ng + “ng + =F (nq MAT+HAT) . (3.1) 


We follow the usual procedure. 


nq nq 
implies 1 
_ 1 (2A (AtnB) 72 
= 2 [2AtaenB) 2 (3.2) 


This describes a "valley" in the n,q plane's first quadrant 
for which the value of q is optimal given any value for n. 
Note that when n is equal to 1 the optimal value of q is 0° 
To confirm that the valley is a valley of minimums, not a 


ridge of maximums, we investigate the second derivative. 


2 

aK Ax nB) 

r2 — + — PIO! 3 
g ng ng 
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Now substituting equation (3.2) into equation (3.1), we get 
the cost as a function of n only, with optimal (minimum cost) 

‘ q chosen for each n by equation (3.2). The equation for the 
cost in the "valley" is 


A. 


ee. j2daens) ]? _ awa rs ICKT 
2 ne 2 2 


_ A(A+nB) (IC) Sy 


NIP lH 


[2 (At+nB) ] 


1 1 
(2AC) * (AtnB) 7 = (n-1) AEE ; 


i} 


Now to find the optimal value of n we consider the derivative 





a 
i. denen | ie 
on 2 : 
2 (A+nB) * 
. Solving for n we get 
eee mee 
4 AICT 


but this is the worst value of n, since 


zt 
2 
92K -(2X1C) 2B 
ary Ts ak 
on oy 
4 (A+nB) 


implies we have a maximum. That is the valley has a high point 


and decreases toward the right and left from the point where 





arct? 


This value may be positive or negative and the portion of the 
"valley" we are interested in may slope downward to the right 
or left or to both right and left. Figure 5 shows a contour 

map of a typical K over the first quadrant of the n, q plane. 


The unconstrained minimum is at n = ~ where the "valley" is 


eye) 








low. Point P is the highest point in the 
"valley". 


FIGURE 5 


CONTOUR MAP OF K OVER n, q PLANE 


Note: The "valley". is between two arrows marked 
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asymtotic to the q = 0 line. This is consistent with the re- 
sults of the initial phased delivery model. of section 2.2. 
Now we introduce the following constraints 


n 


lv 


x 


q m, 


lv 


and 

w<eng=-Qs<w. 
W, w and m are chosen as discussed im the previous. section 
and may represent the most strict of several verbalized con- 
straints. Figure 6 shows a typical valley and constraint 
set.. Feasible solutions must lie with the "crescent" area 
and the valley is the line of optimal (minimal cost) values 
of q for values of n. The peak in the valley is not shown, 
since it will not be necessary to. know its: exact location 
which may be to the right or left or withim the constraint 
set. 

In. order to fully describe the surface of K values over 
the n, q plane we need to consider three additional things. 
First, how does K behave on a constraint boundary with 0 = 
nq fixed? To find out we substitute @ for nq in equation 
(3.1): ta get 


x = AL, mBA , ICQ _ nar , ICAt 


Q Q e a —_- 
The only variable terms are a and SEAT 
Therefore, K is monotonically increasing or decreasing with 
respect to n dependent upon whether 5 — 554 is positive or 


negative. The second consideration is: the behavior of K over 


a constraint boundary with q =m. Substitute the constant m 
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Pa tal 


for q in equation (3.1) to get 


x= AA, BA, PCnm _ Wendt, Wert 
~ Tm ie ae Pe HX, 


For feasibility m2 At. If m= Aq then K decreases monoton- 
ically as n increases. If m > At the terms of K which depend 


on n are of the form Cc, /n + Co n where Cy and C, are positive 


constants .. 
_ AA 
a 
and 


Cy = = {m-AT) 


This implies that K has a minimal value for some n such that 
O<n* <~, Specifically 
2 
—— [ 2AA | 
i TGn (m-AT) S| if 
Unfortunately, this point is not necessarily in the “valley". 


Finally, a quick look at equation (3.2) im the form 
i 


re ta < pA taae) . 


shows that the "valley" has higher values of Q for higher 
values of n, since n is always positive. This insures that 
the "valley" rises, relative to lines of constant Q value, 
such as the non-linear boundaries formed by the 0 con- 
straints. 

We now have sufficient facts to determine the optimal 


solution. Let us marshall the facts in order to see clearly 


where the optimal values lie. 





First, for any given value of n, the valley contains 
the unconstrained optimum value of q. The value of K as a 
function of Q only in monotonically increasing as we move 
away from thé valley. Therefore, for any value of n for 
which the valley is outside the feasible solution set, the 
optimal gq is on the boundary "nearest" the valley. 

Second, the valley slopes downward toward at least one 
end, and therefore the section of valley within the feasible 
solution set has its minimum at a boundary. 

Third, the first two arguments insure that the optimal 
values of q and n lie’ on a pandas constraint. 

Fourth): since we: know that along constraints of con- 
stant Q value K monotonically increases to ¢€ither the right 
or left, the opt imun: cannot be on such a constraint. 

A 

Now we can cohol use tbat the optimum feasible solution 
je on the constantsas> 1 with:w sg < Wand.q@2 m, or on 
the constraint gi=th with w/myn x W/m, andn > 1. 

We will need “to compute and compare the cost of two 
solutions to find thé optimal. policy. First we consider the 


candidate. optiimumaiith n’ =: d,c-and’ 1 





eSVLFi me a ’ 2 
2, (A+B) 
a Goi L SECs: ok a ; 
If q <~maxdw,m) ityen }!o=:max'(w,m). IF gq > W then q' = W 
otherwise q' = q. Now compute K" for n' and q'. Then . 
eompute the ober: candidate optimum with q" = m; and 
Sil oJ om E : 5 ' é 


n= Lreateeneyds. 2 





ee ee 


TF n < max (l,w/m) then n" = max (1,w/m). TIE n > W/m then 


n" = W/m. Otherwise n" = n. Now gompute K" for n" and q”. 

The minimum of K' and K" is K* and n* and q* are the 
corresponding n and q values. 

It is interesting to note that, in general, the opti- 
mal policy is either the conventional gingle increment pol- 
icy, or else the opposite extreme; order the minimum allowed 
increment size and the maximum allowed total quantity. How- 
ever, if the Q constraints are severe or m ig comparable in 
magnitude to i} and t is very small, this comment may not 
apply. ) 

3.4 Full Control of Delivery Size and Interval 

Auume we have full control over n, Tt, and gq. Based on 
prior models, we anticipate the need for a constraint on 
maximum order size. That is, we will require nq to be less 
than or equal to some maximum order size M, Since this 
model is deterministic,we will have each cycle identical. 


Now order cost per unit time is 


AL, BAL 
ng q 


and holding cost per unit time is 


= (ng-nAt+At), 


so that 


— a (nq-—naAt+AT) .. 





Let us first consider the optimal value of Tt. 


a ad Fris~ “H+ 

for ‘all n > iL. 
This implies K is linear in t and inspection shows that 
’ ‘4 ; ag rere 

K decreases as T increases. However, for feasibility 


q 2 \T so that we can determine t* in terms of q*. That 
ii Br aiieys. sams, ~ 

is, we have reduced the model to two variables, and 

T* = q*/X. This allows us to substitute q/A for t and K 

becomes 

~~ : veewmen 

-AA, BA, IC. 

nq q 2 


Now we see-‘that since K decreases as nq increases and nq <M 
n*q*. will bé equal to-M; and we can determine ‘n* in terms 4 
6f q*..-We can now substitute M/q for n and K becomes 


8 WC stibeten |Looe ay -BASSe q 
eee ee ee 


** = 


The first term is a constant,so the variable cost K which 


we wish to minimize is 


and we see that we have reduced the model to a single de- 
Cision variable, q. This is immediately recognizable as 


being of the simple lot size form with 


2BA 
Se 
Interpretation .of .the optimal policy is now simple and 
straight foward. The optimal lot size is computed, using > 


the Wilson simple lot size formula with only the B cost 


term used as "order cost". A is spent at intervals of M/X 
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and its cost contribution is determined by the external 
policy decision never to order more than M. ‘The optimal 
interval between deliveries is determined by the relation- 
ship tT* = q*/\. The optimal total order is M and therefore 





CHAPTER 4 


CONSTRAINED, PHASED DELIVERY MODELS WITH | 7 

; UNEQUAL DELIVERY LOTS AND UNEQUAL INTERVALS 
4.1 enatal Dideucaion 

The main application of unequal delivery snevement or 
interval assumptions is to model single period non-steady 
state systems. Such delivery, obviously, could never be 
truly optimal in a steady state. However, we may want to 
accept unequal delivery increments to minimize lead time 
when producer has an increasing production rate as a re- 
sult of having a partial quantity on hand, or applying 
additional resources to the item's production as other con- 
tracts are phased out. While the following models presume 
such an increasing rate, it would be relatively easy to mod- 
el a decreasing rate if desired. The first model makes no 
assumption as to the relative size of delivery increments. 

These models have limited applicability and the more 
basic question of variable production rates is taken up 
in section 7.1. The first and most general of these mod- 
els is considered in some detail and shortages are allowed. 
The latter two models are formulated and a solution tech- 
nique is given but are not fully explored since the assump- 
tions are exceedingly restrictive. 

It is important to recall that in steady state systems 
the optimal procedures result in repetitive identical cycles, 
and any forced deviation from identical replication will be 


nonoptimal.. That is, given completely free choice, one 
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we 


pt Let, 


would never voluntarily choose unequal delivery increments 
or intervals except as a response to prebabilistic influ- 
ences or changing parameters. This means that lewalixaallilyy 
this set of models should only be applied when outside con- 
ditions require their use. 
4.2 Fractional Delivery at Fixed Unequal Intervals 

This series of models is an outgrowth and extension 
of an exercise in Hadley and Whiten at Demand is deter- 
ministic. It is initially assumed that no shortages are 
allowed and that delivery will be made in two increments 
qy and Apr which may not be equal in size. More specif- 
ically, a fraction fi of the total quantity ordered, 0, 
will be delivered first, and then after an interval T the 
remaining fraction fo =l]- f, will be received. Assume 
temporarily that 1t ‘fy the Wilson model optimal cycle 
time. Now q, = £0, qd, = £42 and gy + dy = Q. Since n= 2, 
A+ nB is a constant and A' = A+ 2B -widl be used in this 
model to represent this total order cost. After the basic 
model is analyzed, a simple shortage model will be consid- 
ered for the two pant delivery conditign, and then the ex- 
tension of the no shortage model to n delivery segments of 
size Gps Ign eer G will be discussed. © and b are the 
decision variables. It appears necessary to investigate 


three cases: 


q, = £0 > Mt > £, > r~vv see Figure 7A 
q, = £0 = a1 S> f= t/t see Figure 7B 
qd, = £0 « At >, < t/T see Figure 7C 





1. Chapter 2, problem 2~65. 
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FIGURE 7A 


TWO INCREMENT MODEL 
7 ' G, = £)Q > At 





QUANTITY 








FIGURE 7B 


TWO INCREMENT MODEL 
q, = £0 = AT 
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FIGURE 7C 


TWO INCREMENT MODEL 
q, = £,0 < AT 
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However, it will not be necessary to distinguish between 
these cases during the first few steps of the analysis, 

but reference to all three figures will aid in following 
the analysis. Holding cost per cycle is proportional to 


the area under the on-hand curve. Consider the area under 


the on hand curve during the time t. It is a trapezoid 


with an area of 
is + is aad 
tt ] 
2 e 





Similarly, the area below the curve between the time oN 
is delivered and the end of the cycle is a trapezoid whose 
area is given by 


(T-+1) | See) 2) | , 


Holding cost per cycle is IC times the sum of these two 


areas, and equals 


TE [elabe2e,o-at) + (Pet) (2bHA (t-1))] 


or 


= (2£,ATT+2_T+OT-2ATT) . 


Dividing by T to get holding cost per unit time vields 


= (2£,AT+2b+Q-2AT) , 
= 252 + IC(£,At+b-AT). 


Now since order cost per unit time is A'\/QO, 


- Ar + 102 , IC(£,At+b-At). (4.1) 


K 7) 5 


Taking the derivatives of K with respect to b and Q, the 
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decision variables, setting the derivatives equal to zero 


and solving yields, 


This implies K is a linear function of b and inspection of 
K shows K is minimal when b is minimal. Now, 


~A'h Ic 
2 2 





implies the optimal value of Q is Qn This is a somewhat 


surprising result since one would probably expect Q* > oe 


However, to continue, it is necessary to consider the var- 


ious cases to determine b*. If fy 2 T/T then b can equal 


0 and b* 0. Now 


ECO 
Ww 
2 


i] 
K* = ee 


QW 





+ TCAt(£,-1), 


but £, - l= -£ so that 


x = a 
K = ICAtE, s Ky 


However, if fy < t/T, there will be shortages unless 
b +£,0 2 At. This implies that b > At - £,Q and the mini- 
mal allowable b is the optimal b* = it - f oe Now the 
optimal values b* and Q* can be substituted into the 
equation for K to get kK*; 
K* = A’) is TCQ,, 
ea os 5 
# 





+ IC(£,At+At-£ O77 AT) 


au 


KO + ICF,A(t-T) 


But, since we assume Ty > t and I, C, fy and X are all 
greater than zero, K* can be written so as to make its 


relationship to Kk obvious; 
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x = ce mi 
K K Icf,i(T,, Th = KY ‘ 


It is now clear that, although Q* = Qe the split delivery 
results in a reduced cost. Under these assumptions, the 
optimal operating policy is the same as for the iinpte 
Wilson model with the addition that b* = 0 or b* = dt - £12 
depending on the case occuring. 

Actually, only one case is required, as can be readily 
seen by reference to Figure 8. Since Q* = % and T* = Th 
the values Gy Age T and Ty - t are fixed and the cycle 


repeats identically. We can now see that we can define a 


cycle for computation purposes to be Pi to P. or equally 


well from P, to Py. ‘That is, we can choose f' = f. if 
and only if fy 2 t/t), and f! = fi otherwise. We then de- 
fine tif tor 1? = Ty - T respectively. Now 
Q* = Qn 
| 0 if f) > t/T, 
= ’ 
nn £12, if fl < T/T, 
and 
* = = t ' 
K Ky IE * AT" os K . 


It is considerably less obvious how to apply the above 
formulas when Tt > Th Figures 9A and 9B show cycles 


graphically for Ty “Gia aT and aq =". & 3T . More 


than one order will be outstanding at times. Define 


C => = nT such that Ty >t" > 0. ft is clear that 


Le) 
+ 
it 


7 is still optimal. It is also readily seen that 


can be used in place of t to determine f' and t' and 
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that K* = K, - ICf'At' will be the proper minimal cost. It 
remains to define b*. It is obvious that the same technique 
of using t" in lieu of Tt gives the pitipneat optimal values for 
D~.. , 

Suppose now, that shortages are allowed and there is a 
cost 1 associated with the occurrence of each shortage and 
a shortage cost rate 7 per shortage per unit time. In the 
lost sales case, T must be equal to zero since it is mean- 
ingless to discuss the duration of a lost sale. ‘The cost 1 
may be interpreted in this case as the profit lost as a re- 
sult of not making the sale. Both 7 and T may be applicable 
to the backorder case. The term 7 is interpreted here as 
the cost of processing a backorder which is independent of 
the time for which the backorder is outstanding. The a term 
is the portion of the backorder cost which is dependent upon 
the duration of the backorder. 

Some preliminary calculation can be made on the effect 
of shortages on holding cost and of the direct shortage costs 
prior to modifying the formula for total cost, K. Define s 
as the maximum number of shortages occurring at one time. By 
the previous arguments reducing the problem to one case, we 
will assume that fy > t/T. There are two possibilities for 
shortages. They may occur prior to one or both deliveries. 
Figure 10A illustrates the case where shortages occur only 
once during a cycle, at the end, since fy > t/T. The number 


of shortages per cycle is s and the average number of short- 


ages per unit time is - s (s/\). Therefore shortage cost 





per unit time is: 


Holding costs sae unit time must be reduced by 
0 (s-$5 ). 

Figure 10B illustrates the case where shortages occur 
prior to both deliveries. The number of shortages per cycle 
is 2s - (q,-AT). As expected, if q, = At, the number of 
shortages per cycle is 2s. This term is multiplied by 7 and 
divided by T, to get the shortage cost factor to add to K. 
In a similar fashion, the applicable cost factor involving t 


is found to be 


ee [s? + (s-(£,0-A1))? | 


The holding cost adjustment is: ? 


5? + (s-(f perca bigee 
eo 
20 


Now there are two cases to be considered. 


Case I- q 2 At and shortages only at end of cycle. 


Case II qi ae: and shortages before each delivery. 
These cases can be solved by adding the appropriate factors 
above to equation (4.1). It is also convenient to substitute 
£,2 for qy to make obvious the prover functional dependence 
upon Q. 
It is no longer necessary to include b in the equations : 
as a decision variable, since s serves the purpose for calcu- 


lation and optimization and the value of b* = -s*. The 
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TWO INCREMENT MODEL WITH SHORTAGES 
AT END OF CYCLE ONLY 


QUANTITY 





FIGURE 10B 


TWO INCREMENT MODEL WITH SHORTAGES 
BEFORE EACH DELIVERY 





equations for K are given below. 


Case I 
Ad . Ico 5? ms . TS" : 

K = ns + = + IC(£,At-st+55) + rae + “20 . 
Case II " 


s*+(s-(£,0-A1))* 1 


ai 20 


(+ FRQs- (0-40) + Fy [8° + (s-(£0-20)7 ] 


To get -explicit results for case I, it is necessary to take 
the derivatives of K with respect to Q and s, set them equal 


to zero, and solve. 


implies 
2,* 
eas) +s cat) (4.2) 


implies 


or 


_ Amts (IC+T) 
oes eee (4.3) 


Calculation of the second derivatives shows that we have a 
minimum. Combining equations (4.2) and (4.3) yields the 


following quadratic in s; 


(q°.% Gidye” + Qenie = Dare. ~ Gay? 


Note that if 7 = 0 then either there is no solution for 


2 


0 <s < ~ or if (Am)” = 2\AIC any value of s is a solution. 
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In words, this implies that if 7 = 0 and v # Q, then either 
there should be no shortages or infinitely many. That is, 
the system should incur no shortages or else go out of busi- 
ness. The lost sales case is characterized by just these 


conditions 7 = 0 and - HO RE Tt # 0 “i. 


_ |Paa(n+re)— (a0) - (At) 


Q* 
TLC 
a 
oi [2AAZS _ Ie(an)* P TA 
(ice) a = (reer “F (rt+T) 
and 
be ae 


These results are identical to those for the single increment 
version, however, the cost K* will be less by ICT, \T. Ex~ 
plicit results for case II could be gotten by the ‘same 
straightforward procedure but with increased algebraic dif- 
ficulty. 

The basic two delivery, no shortage model, is readily 
extended to a more general n delivery, no shortage model. 
Figure 11 illustrates such a model. Modified notation is 


required. Assume n fractional deliveries fy ,fo7--+ +t, 


n 
where , & = ll. When £,Q = dy and there are n deliveries 
q; = Q. T, is defined for 


i=1,2,...,n-1 to be the time between the ith and i+lst de- 
n=-1 

livery, ) 1, < T for a feasible solution to exist. It is 
i=1 


obvious by comparison with the basic model, that 9* = Ls A 
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It is only a bit more difficult to conclude from the defini- 
tion of b and the geometry of the model that 

i 
b* = MAXIMUM (0,ary-ay, A(T +T5)— lay tag) reves fron 


1 


n n 
r = a 2 x -9) « 
1 5 dy aa y 2 


j= 
Similarly an expression for K* can be found. 


} a ' 
a ak - CT ee = 7 (4, 5) . 
= ig fn 27) 


No attempt is made to consider the n increment model 
with shortages allowed here, for several reasons. First, 
the formulation involves too many different cases to be of 
much interest unless severe and unrealistic aceite are 
made. In fact, there are 2n-1l different possibilities for 
distinct cases, since the on hand curve has 2n distinct 
vertices and the zero line may fall optimally between any 
pair of them. The solution procedure would, of necessity, 
follow the technique of assuming a particular case and deter- 
mining if the optimal values are consistent with that assump- 
tion. After all consistent solutions are found, it is then 
necessary to compute and compare the cost K resulting from 
each. The minimum K is K* and the associated optimal values 
are the optimal policy. This is a tedious technique, but is 
usually susceptible to computer solution for any specific 
application. The second reason for not pursuing this further 
is that it is doubtful if the model would have any sigqnif- 


icant applications for a value of n greater than two. It 
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seems far more likely that if the number of deliveries is 
large, their size and the time between deliveries would be 
uniform, or at least subject to small random variation. In 
this event, other models presented provide a close approxima- 
tion and are much easier to use. 

One particular delivery pattern that might warrant fur- 
ther analysis at a later date is the possibility that either 
the first or last of a large number of shipments might be an 
odd size, while all of the other lots are the same size. 

This could occur, for example, as the result of the source 
having a partial quantity on hand for immediate delivery and 
producing the remainder and shipping in equal increments. 
‘Tt appears probable that significant simplifications could be 
made. | 
4.3. Increasing Delivery Lot Size at Equal Intervals 

| It is assumed in this model that a manufacturer with an 
increasing Srodudbien rate has agreed to make shipments at 
fixed intervals, of the entire production since the eae shiv- 
ment. Let t be the fixed interval and Py the production in 


the 1th period. Ds 


7, Ps if and only if i > j Since we as- 


sume an increasing production rate. The usage during each 
period is At and the only parameter over which the manager 
has control is n. No shortages are allowed. See Figure 12; 
Now there exists an i such that pete AT, 


and there exists a j such that ; Dy > jAt, 
i=1 


therefore the system has a feasible solution. 
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FIGURE 12 


INCREASING PRODUCTION RATE: 
DELIVERY AT FIXED INTERVALS 





n 
Define Q(n) = ) p.. 
i=l * 


The on hand level required at the time of the first delivery 


in a series to prevent shortages is - 


a. 
b = max (iAT - ) Du). x 
i j=1 J 


‘ . ; Ad 
Now the order cost per unit time is Oty ! 


and holding cost per unit time is 





n 
iG [ 2)n +b- ) (i-1) tp, |. 
: i 
1=2 
The formulation of the model is then to minimize 


n 
Kin) = stay + re [2m +b-t }) Ip = >, | . 
i=2 





subject to n a positive integer and ) D, 2 nat. Specific 
i=l 

solutions are fairly easy to find when the production quan- 
tities have a convenient algebraic relationship. As a 
reasonable example, assume the production rate is a linear 
increasing function of time. That is, assume ¥ = a+ a't 
where a' is greater than zero and a+ a' is greater than 
zero. For convenience,we will use t as our time interval 
so that pi = at ia'. Figure 12 shows an example of the 

on hand value resulting from this delivery pattern. The 


decision variable is n, the only parameter under our con-~ 


trol. The total order quantity Q is a function of n . 


denoted Q(n) where 





BR fens 


ata'+a+t2a' +... ¢+ atia' +... + a+ na’ 


i] 


Q(n) 


il 
BD 
+ 
=. 

~~ 

P 

l 
B 
ae 


To prevent shortages, a quantity b must be on hand at the 
time of first delivery. The on hand decreases at each 
successive low point, until a+ ia' > At. Lety=i-l. 


Then y is an integer and 


b= vty - ya - S¥jyet) 


It is necessary that the order quantity Q(n) be greater or 
equal to the demand during the time over which it is 
delivered. That is, 

QOtm) 2 nar, 


or 


a + a' (n+l) » ie 
2? e 


which can be rewritten 


2AtT - a-a' 


n > = 


We know the order cost per cycle is A, so that the order 


cost per unit time is 


iN} 2Ad 
Q(h) ni@atanra’) ~° 


The holding cost per cycle can be seen from Figure 12 using 


previously described techniques to be 
A [ doin) z= (<(at2a") + 2t{at3as) + ane + DTUAtR FE | ues 


+ (n-1) t(ata'+(n-1)a") ) + bt | " 
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so that holding cost per unit time is 
n-l n-1 
Ic 7 ICT?) ee 
= Q(n) + ICb ae i(ata') + ) Par 
1=1 11 
ICb is independent of n and can be dropped. We can now com- 
bine these terms to get total variable cost 


2Ar TCan ICa’ (n“+n) 


iS n(2ata'nta') = 2 a 





ICAT Za'n*+3an¢6a!n~Jaq2a! 7 
3 2atna'+a' 


Setting the derivative of K with respect to n equal to 

zero, yields a 5th degree polynomial with exactly one sign 
change which by Descarte's Rule implies exactly one positive 
root. The second derivative is positive, so there exists 

a unique n > 0 such that K is a minimum. There is no 
algebraic exact method for computing the root of a 5th 
degree polynomial, so the Newton-Raphson method or some 
other convergent iteration technique must be used to find 
this value of n. Call this value n'. If 


; 2\T-a-a' 
2 ' 
= a 


n 
it is feasible and therefore the desired solution is n* = n',. 
Lf 


2r\T-a-a' 
Se 


then n* is the smallest integer larger than or equal to 


2\T-a-a' 
at a 


since K(n) increases monotonically with n for all n> n'. 
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4.4 Equal Delivery Lots at Decreasing Intervals 

In this section, it is assumed that a producer has an 
increasing production rate and will make a shipment of a 
fixed size at unequal intervals. For clarification, an 
example is given. United Widget makes gismos for our com- 
pany. When we contact them, we are informed that they have 
a government contract in production now, but that they can 
produce 100 per month for us now and will add a 100 per 
semi production line each month for 8 months as the govern- 
ment contract expires. Gismos are much cheaper to ship in 
carload lots of 200, so he offers to ship each carload lot 
of 200 as soon as it is completed. We would expect ship- 
ments at approximately 1.5 months, 2.3 months, 3 months, 
3.5 months, etc. 

The general problem can be formulated as follows. We 
get a quantity q at intervals Ts where the increasing rate 
implies that Ts > T3 if and only if i < j. Our decision 
variable is n, and the total order quantity will be nq. 

n-1 
For feasibility nq > 3h where t,) is the delay until 
the first shipment is received. Figure 13 shows an exam- 


ple of a cycle of this model. It is obvious using the, by 


now familiar, techniques that 


n-l 
K(n) = AA 4 2089 _ icq] J it : 
nq 2 j=l n= 


For convenience,we will use difference equations rather 


than taking derivatives. That is, define 
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AK = K(n+l1) - K(n). Our optimal integer value of n is the 
smallest feasible value of n such that AK is positive. 


Find the smallest n such that 


n 
3 ~Ad Ic CK 
= ———— + + - = .) ct. & = 
q (n° +n) i= 


This can be rewritten as 





2 A # 2A) 
(n“+n) | 1 -— } “| > ‘ 
nq 11 q*ICc 


It is clear that the left side of the inequality is a mono- 
tenically increasing function of n since the T; are de- 


creasing. Now we call the smallest integer value of n that 


satisfies this inequality n'. For feasibility we must have 
: n-1 
r 
me= } t, w 
1 4a=p * 
so that if 
n'=1 
rv 
n' 2+ ) i 
2 i=o0 * 


5 
A 
Q{_> 
a 


i=0 


then n* is the smallest integer value of n such that 


" n-1 
ny ee i= 
4 i=o + 
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CHAPTER 5 


SPECIAL CONSIDERATIONS IN 
PHASED DELIVERY MODELS 
5.1 Shortages 

The problem of incorporating shortages and shortage 
costs into incremental delivery models can be exceedingly 
complex or, in special cases, quite simple. There are 
several techniques for approximation which may be nearly 
exact under certain conditions. In this section, a tech- 
nique for exact solution of deterministic models with short- 
ages will be described and three approximation techniques 
will be presented. Reference to earlier models and simple 
examples will be given as appropriate. 

Figure 14 shows a typical graph of on hand level de- 
livery for a cycle of n incremental deliveries, resulting 
from a single contract. It is obvious that the on hand 
level can be lowered or raised in relation to the zero 
stock level to balance holding and shortage costs with 
each other and with order costs to achieve a minimum total 
cost. The problem in formulating an exact model is based 
on the fact that there are 2n - 1 different forms for the 
shortage cost, depending upon where the zero line cuts the 
on hand curve. There are 2n points connected by the on hand 
curve and the zero line may lie above 1, 2, .... 2n - 2, or 
2n - 1 of them. Call these levels 1, ... , 2n - 2 and 2n - l. 
It will never be optimal for the on hand curve to lie en- 


tirely above or below the zero line if infinite and zero 
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shortage costs are not allowed. In these cases, the on hand 
curve would "sit on" the zero line or no inventory should 
be kept respectively. The level at which the zero line 
cuts the on hand curve determines a value of s, the number 
of shortages. The functional form of holding cost and 
shortage cost as a function of the number of shortages per 
cycle is different for each level of the zero line. There 
“is a one to one correspondence between these levels of the 
zero line and ranges of shortages which may result, so that 
for each functional form there is a feasible range of values 
for s, and for each value of s there is a functional form. 
This allows us to define K(i) as the minimal variable cost 
per unit time for the ith range of s where i=l, 2... , 
2n - 1. K(s) as a function of the number of shortages is 
then defined by K(s) = K(i) for all s contained in the ith 
range. Figure 15 shows such a cost function. The problem 
«ig remarkably similar to the incremental price discount 
models found in the literature. The computational procedure 
is virtually identical. First, compute the optimal policy 
Ps for Gach K(i) idm 1, a.., 2H - 1. The policy Py in- 
cludes a value for shortages which may or may not lie with- 
in the allowable range for s. For each of the policies Py 
for which S5 is included in the ith range for s,compute 
K(i,p;)- The optimal policy p* is then the Ps for which 
K(i,p;) is minimal. 
Less it sound overly simple to solve such a problem, 


recall that the 2n - 1 functional forms may be very complex 
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and disparate if the Ts and q, are unequal. The zero line 
may cut the on hand curve at 1, 3, 5, 7... different points 
and the problem of specifically defining the holding cost 
and shortage cost functions may be very difficult. The 
procedure is straightforward, but may be very difficult in 
application without the aid of high speed computers. 

The first approximation is relatively simple minded 
and applies exactly to models which result in q, = 4 and 
T; = Tt for all i and for which q = At. We assume a model 
comparable to that in section 3.4 as shown in Figure 16. 
Note that the same number of shortages will occur at the 
end of each period t. We see this approximation is good 
whenever gq-At and n are small. If we let q-At = e,then 


the approximation overstates the number of shortages by 


nine and results in an "optimal" policy which has a 


value for s* which is too small, so that our policy will | 
be cingekveriaes in the sense of noiiiewing a hide level 
of service than forecast. It is interesting to calculate 
specific results for the shortage version of the model in 
section 3.4. Recall that gq = At will be optimal. 

Assume the source of an item agrees to produce at a 
rate equal to \, the demand rate. The values of A, B, C, 
I and A are known. Shortages will be backordered and the 


cost per unit per unit time of a backorder is 7. The cost 


per backorder independent of time is mt. Our decision vari- 


ables are q, n and sj; the number of shortages prior to 
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SHORTAGE MODEL 


APPROXIMATION ONE 
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receipt of a lot. Now order cost per unit time is 


holding cost per unit time is 
Ic (q-s)? 
2 q : 


and shortage cost per unit time is 
2 





as TAS : 
2q q 
Now we see that 
52 oD 
Ad Br oe T Ss TAS 
K=—+—+tH+ 2 c + +—, 
ng q ial a q 


It is immediately obvious that K is minimal with respect to 
n when n is unpfdaee. Since no real system would allow an 
infinite order, we will assume nq <M is imposed. For ex- 
ample, top level management's policy may be to never order 


more than two years worth of stock. Now the variable cost is 


32 
Kos LIS “- ie - TCs + (IC+n) on + ms P 


The optimal value of n is 


Using the usual technique, the derivatives with respect to q 


and s yield 


2Bh (ic+n)s | pue z 
a. So = 


= gic-m 
IC+t 


a mic ca » 2 
respectively. The second equation can be rewritten as 


= (Ic+T)s | 
q Tc 
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It is interesting to note that when 1 = 0 the ratio of ‘the 


time stock is on hand to the time stock is out ‘equals 


fg-s)A _ 
Si 7 old 


the ratio of the shortage and holding cost rates. Setting 
the two equations for q in terms of s equal and rearranging 
we get 


(a7 


+ TIC)s*-+ 2nms + (mA)? ~ 2B\IC = 0. 
It is obvious that if 7 = 0 and 
(ma)? 2BAIC , 
then s = 0 or ~, and it is optimal to either have no ‘short- 


ages or to go out of business. So we require 1 > 0 and solve 


for the optimal values of q, s, and n. We find 


1 
ar = [(ESE) (BR - = ll 


IC (T+IC) 
gh nx q* TCA 
(T+IC) 
and 
n* = M/q* . 


These results are very: closely related to the single lot 
model's results. The optimal values for s and q are iden- 
tical, except that the lot order cost, B, is used as an order 
instead of A. As we have seen before, we can use the stand- 
ard formulas and then order a sufficient number of increments 
of the prescribed size to make a total order of M. 

The second approximation technique assumes that the cost 
of shortages is high, relative to the other costs of the mod- 


el, particularly the holding cost, so that shortages will be 
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rare. Shortages are assumed to occur only prior to the first 
delivery in a cycle. Specifically, it is assumed that 

q, 2 AT, + s and the approximation will be exact if, in fact, 
this condition holds for the optimal values. That is,if 

q* > At* + s*. This approximation cannot be applied success- 
fully when we have full control over delivery conditions since 
then q* = At*. Reference to Figure 17 will assist in recog- 
nizing that this approximation results in under stating short- 
ages in the model. It is not a conservative policy, since 

the actual number of shortages and duration of shortages will 
exceed the model forecast and result in reduced service. 

For example, a shortage model similar to the model of section 


3.2 where q and t are fixed would be approximated by 


we 
K = a ae ZC (ng-ndtT+\T-28 a 12,4. : 
nq q 2 nq nq 2ng ‘ 


The only decision variables are n and s since q and 7T are 


fixed. 

2 A 
oK _ -Ad Ic Cs _ 7™srX _ TS 
oA oe ee ae ee 

nq nq nq 2n'q 
and 
8K _ era + EOS, MT 2 = 2 
s nq nq nq 
imply 
s* ao asl : 
IC+T 

and 
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SHORTAGE MODEL 


APPROXIMATION TWO 
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The final approximation technique to be presented here 
is designed to fit systems where q > At and both q and T 
are small while n is large. It is basically an approximation 
which reduces the "saw tooth" form of a rising on hand curve 
to a straight line,as in the continuous finite production 
rate model. Figure 18 shows a typical cycle of incremental 
deliveries with several alternate approximating straight 
lines superimposed as broken lines. Our problem is to choose 
one of the infinite number of straight line approximations. 
Line A will result in overstating holding cost and under- 
stating shortage cost. Line C will understate holding cost 
and overstate the time dependent aspect of shortage cost. 
All understate the number of shortages. The policies result- 
ing from the use of a straight line become more conservative, 
in the sense defined earlier, as the approximating line moves 
fue A to B to Cc. In the limit, as q and T go to zero while 
n goes to infinity, a mid-point approximation such as line B 
is asymptotically exact. It is not unbiased, since short- 
ages are always understated. On the other hand, line C has 
the advantage of being the most conservative approach and 
also provides the "best" forecast of shortage costs, since 
it understates the multiplier of 7, while overstating the mul- 
tiplier of tm: We can wait until after the’ calculations are 
made to decide which approximation to use. 

Figure 19 shows a single cycle of the approximating 
model. The slope of the rising on hand curve is g/t - i. 
Recall that if either q or T are variable, then q = AT and 


this model will not apply, so our decision variables are n 
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FIGURE 19 


APPROXIMATE SHORTAGE MODEL 
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and s- Also recall R = q/t. Now in the approximating model 
consider the areas marked Al, A2, A3, and A4 on Figure 19. 
The holding cost per unit time is 


(A1l+A2) 


= TS. 1 


The shortage cost per unit time is 


(A3+A4) * TS 
Ce. ae 


Therefore the applicable variable cost per unit time is 


¢ gS Ftc 


: ee 
R-} 27 R Ts“/ R ue 
q*q ” ng na( R ) s| =x) * Ing =x) ° na * 
Using the usual process we can determine n' and s', the op- 


timal values for the approximating model, which are 





cl 
Law x [Ss La) at. (wa)? |= 
Thom. ee ae 


and 
v= Cet) (Ca). 
IC+n 

Now we must decide which approximation to use. n* = n' 
for all three cases but in setting the reorder point we can 
use s* = s' + a where a is an adjustment to yield a more or 
less conservative policy. Increasing s* increases the re- 
order point and reduces the actual occurrence of shortages. 
a = 0 is equivalent to using the approximating line A of 
Figure 19 while letting a = At is equivalent to the more 
conservative policy using the approximating line C of 


Figure 19. Any value of a can be used. The range 
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0 < a s$ AT Corresponds to the more reasonable set of line 
approximations between line A and line C and the asymptot- 
ically exact value for a is \T/2. 
5.2 Probabilistic Parameters 

The introduction of a probabilistic parameter intro- 
duces considerable complexity into any inventory model. 
This is particularly true of models which have realistic 
features, such as shortage costs, and are designed to be 
applicable to real inventory systems. The introduction of 
a demand distribution into a phased delivery model may be 
exceedingly difficult unless simplifying approximations are 
made. It is assumed that the reader who has gotten this 
far into this paper, is well acquainted with the basic 
single delivery probabilistic models for continuous review 
and periodic review. For the reader who is not, it will 
probably be necessary to study these basic probabilistic 
models in an inventory text such as Arrow, Karlin and 
Scarf [2] or Hadley and Whitin [1]. The following will be 
basically a comparative discussion and consider the modifi- 
cations required rather than derive probabilistic models. 

The modifications required are fairly easily divided 
into three groups. Incremental deliveries effect lead 
time, on hand levels and the frequency with which a risk 
of shortage is taken. Essentially,the desireable effects 
are based on the reduction of leadtime and stockage levels, 


while the deleterious effects arise from a more frequent 
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risk af shartage. The three functions are taken up in 
order of increasing complexity. 

It is quite probable that phased delivery models 
will result in earlier deliveries of the first units or- 
dered. Unless constrained externally to the model, q* 
will always be less than the Q* which would result from a 
single delivery. This will usually mean the delivery quan- 
tity q will be available for shipment a shorter time after 
the order is placed than will be Q*. We have defined Tt, 
as the delay prior to receipt of the first delivery of 
material. Since safety stock is monotonically nondecreasing 


with respect to T,, we can expect a reduced safety level 


0 
(and therefore reduced cost} due to incremental deliveries. | 
If the safety level is large enough to make the possibility 
of shortages after arrival of the first increment remote, 
then simply using To for the lead time will take care of 
this change. If the problem of shortages later in the 

cycle is significant, then it would be difficult to get 
exact results but the distribution of demand over Tor 


+ Ty, T+ 1, +T 


T9 0 1 
and used in the probabilistic model. 


2" and so on, would have to be computed 


Fortunately, the problem of determining the reduced 
level of on hand stock is simplified by the fact that the 
reduction in net inventory due to the incremental deliv- 
eries will be known precisely. The problem is that the re- 
duction in holding cost is based on the reduction in aver- 


age on hand stock and not on average net inventory. Once 
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again, it is clear that if shortages are rare and it is 
extremely improbable that shortages will occur after the 
first lot is received, there is no problem. On hand and 
net inventory levels will: be identical over the period of 
time for which the deliveries are placed. This means we. 
can simply subtract from the probabilistic computation of 
holding cost,the fixed saving due to incremental deliveries. 
If the probability of shortages after the first lot is 
received is significant, it will be necessary to calculate 
the expected value of on hand separately for each time per- 
‘iod. That is, integrate over the time between order and 
first delivery, and between each pair of deliveries sepa- 
. rately, using the final distribution of inventory position 
aun one interval as the starting distribution of inventory 
position for the next interval. This is not particularly 
tractable even when the distribution assumed is simple. 
However,: as soon as it becomes obvious that the probability 
of shortages is insignificant, the technique in the previous 
paragraph can be used for the remaining part of the cycle. 
The increased frequency of shortages is a very important 
consideration when q is only a little larger than ATar Or; 
more generally, when the sum of the quantity received in the 
first few deliveries is not much larger than the demand dur- 
ing that period.. When qy >> AT, and receipts occur more 
rapidly than demand is likely to be generated, shortages will 
probably occur only once per cycle. However, there will al- 


ways be some more risk of shortage in an incremental delivery 
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model than in a single delivery model because the build up 
af assets is phased and not instantaneous. 

When special conditions exist, there may be fairly 
simple ways of handling shortages. For example, suppose 
we have r units on hand at the time an order for n incre- 
ments of size q each, at intervals T; as defined earlier 
is placed. Then the expected number af shortages given 
the distribution of demand over a time period t is f(x,t), 


is given by the following sum of integrals. 


f (x-a®) £(x, ty) ax a (xa) £ (x, Tgtt, Jax aes 


af’=srtl a'=rtl+q | 
| 
00 i i | 
+ Sf (x-a)fi(x, J 1.) dx +... 
a =r+I+iq _ 


In general, it will always be possible to consider the 
problem in n+l time intervals; from the placing of an order 
to delivery 1, delivery 1 to delivery 2, ... , delivery 
n - } to delivery n, and delivery n to next order. The 
first n intervals correspond to T; © = Dy, wes pte Tis 
Define t, as the interval between the last (nth) delivery 
and the next order placement. 

In this model, when lead time is stochastic, only 
the first interval varies in length, since we assume the 
Tt, are fixed for alli=l,...,n- 1. This introduces 
a real problem, since the distribution of lead time de- 


mand is based on both the distribution of lead time and 


of demand. The lead times for all deliveries except the 
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first, have the form of the sum of two random variables; 
one the lead time demand for Tp and the second, the demand 
over a fixed time period. However, we presume that if To 
is large relative to the tT; Wwe can use the lead time de- 
mand distribution for all deliveries simply by replacing 
To with the appropriate sum of Tye If Ty is small rela- 
tive to the T; we can use a simple demand distribution for 
the various fixed intervals using the expected value of T,. 
With this assumption, we can compute the desired hold- 
ing cost and shortage cost terms by summing the results for 
the holding and shortage cost terms for each of the n+ 1 
intervals computed as follows. For the jth interval j = 0, 
1, 2, «.., nN use a constructed reorder point of 
rye lai 
where Ig = r, the conventional reorder point. Use 
i 
i=0 7 
as the parameter for mean lead time. For limits on the 
4th time interval integration use 
j-1 
T. and le 9 
i=0 + i=0 * 


as lower and upper limits respectively where 


Mi 
4 
p- 
Vil 
oO 


This will usually be a tedious and lengthy computation, 


but each integral in the sum will be tractable if the 


92 


corresponding single delivery model would be tractable. 
If not, expected value models would have to serve. 
5.3  Reducible Models 

Whenever we can delay determining the size of the ith 
lot beyond the time the contract is placed, we can reduce 
the model into stages. The total order quantity is deter- 
mined by one model. The optimal shipment lot sizes can be 
determined in stages, using a second mode] at the time the 
decision must be made, say one month prior to the delivery. 
This is a typical contractual situation. Often annual con- 
tracts are placed for an approximate quantity and calls are 
placed each month for the quantity needed. There is a po- 
tential problem in determining the duration of the contract 
but generally this can be ignored since it will be set by 
policy. The real problems are in determining the optimal 
interval between orders and the optimal quantity to be or- 
dered. The optimal interval between placement of call or- 
ders against the contract can be found using a modified 
version of the model used to determine the optimal cycle 
length when operating on a periodic review basis. Having 
determined the optimal interval, the optimal quantity and 
safety level can be calculated using a cqnventional per- 
iodic review model with the incremental lot order cost 
used in place of the contract placement order cost. These 
problems are well covered in the literature. 


In a forthcoming paper sponsored by the U. S. Army 


Material Command, Deemer and Hockstra [3] have developed 





a model of this type in detail. They call it a decounled 
model since the costs can be decoupled into those depen- 
dent upon the contracting cycle and those dependent upon 
the safety level and the two problems separated. They 
assume a fixed delivery cycle of one month. This model is 
useful if the control over each delivery quantity is ex- 
tensive and the decisions can be made as time progresses 
rather than in advance. Because of the paper referenced, 


this problem. will not be further discussed here. 
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CHAPTER 6 


STOCHASTIC DELIVERY QUANTITIES 

Delivery quantity may be stochastic for a fairly ex- 
tensive class of ans, Chanecentaadl by the fact that the 
cost may be signi tloantty lower if the source is allowed 
some latitude in the quantity to be provided. One example 
of this class would be a valve machined from an expensive 
casting, which may have flaws which are first discovered 
during machining. It is frequently necessary for vroduc- 
tive efficiency to cast all of the blanks at one time, in 
advance. There will then be a random number of blanks 
found to be flawed during machining. If the contract is 
for an exact quantity, the problem of how many blanks to 
cast is an interesting problem itself. It is costly to 
have too few perfect blanks and have to cast some eore, 
losing both time and set up costs. It is also costly to 
cast too meaty blanks for disposition as surplus. Since 
these pawts will be passed on to the purchaser, it will 


often be mee: for the inventory manager to negotiate 


a contract to accept some range of quantity (say 400 + 10%) 


at a fixed cost per unit. This will result in procurements 


for which the delivery quantity is stochastic. The deliv- 


ery quantity need not even have an expected value equal to 


the quantity ordered. 
For example, if the cost of producing surplus blanks 


is very small, (or their salvage value is high) relative 
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to the cost of running short of blanks, it will be desire- 
able to produce many blanks and the expected delivery quan- 
tity will be close to the upper end of the range. On the 
other hand, if the blanks are expensive (make up a large 
sreeoreten of the einai cost of the item) and have little 
or no Salvage value, then relatively few blanks will be 
produced and the expected quantity delivered will be in 
the lower part of the allowed range. The reader will note 
that the producer's problem may be treated as a modified 
newsboy problem. | 

Suppose one has determined an optimal operating policy 
based on a model which assumes delivery of the exact quan- 
tity ordered. When placing the contract, the price quoted ‘ 
for an aseaet number of units is $100 each, but the con- 
tractor otters a price of $98 each if you will accept a 
variation of plus semis 10% in the delivery guantity. 
He is willing to do this for reasons similar to those 
mentioned eerie: Is our computed policy optimal, or 
even near optimal, now that we have this additional option? 
To find aie we will investigate the revised optimal Lt 

We know that the cost function K is relatively Ciek te 
the vicinity oa aes minimum, so a change of plus or minus 
108 ‘a oie veetiict ghaneees will probably not increase 
costs auerietentiy to offset the reduced unit cost. The 
ratio of the cost of a nonoptimal Q to the cost of O* is 


given by 


(f+) 
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For a nonoptimal Q, which is off by a factor of 2 in either 
direction, there is an increase of approximately 25% in the 
cost. For a nonoptimal Q, which is within plus or minus 
10% of the optimal value, this ratio is less than 1.005. 
This implies an increase in operating cost of less than + 
of 1% compared to a 2% discount on unit cost. If the an- 
nual operating cost is less than 4 times the value of annu- 
al demand, it will pay to accept the offer of the variable 
quantity at the lower cost per unit. However, there will 
be some increase in operating cost, due to the uncertainty 
about the quantity that will actually be delivered, which 
has not yet been considered. 

Let us divide the discussion into two parts. We will 
consider the increased cost due to uncertainty later. 
First, assume that the lead time demand is always less than 
the smallest possible delivery quantity. Then we will have 
no orders outstanding to contribute uncertainty at the time 
an order is placed, so that the reorder point computation 
can be made as if there were no variation in Q. However, 
by adjusting the quantity, we may be able to further de- 
crease costs. Assume that the density function £(Q) for 
the delivery quantity Q based on an order of size Q; is 
known. We are actually supplying one of the parameters of 
the distribution £(Q) when we determine Q'. Now, we can 
choose Q' so that the expected value of K is minimal. Let 
r be the smallest allowed ratio of Q to Q' and r' the 


largest. We want to find Q' such that 
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r'qQ'! 
f K(Q)£(Q)dQ 
rQ' 
is a minimum. 
Suppose, for example, that . 
K(Q) = a/Q + bQ, (roughly the form of Ko 
and that Q is distributed uniformly from .90' to 1.1 Q' 
(+10%). We want to choose Q' so that the expected value 


of K will be a minimum. That is, minimize the following 


integral with respect to Q': 


1210" 


1 
K' = f £4 bQ eer do 
-9Q' Q -2Q 
v a! ' 1 
K(Q') = or +t b'Q 


where a' = 5a 1ln(11/9) = 1.0024 a, and b' = b. The optimal 
value for Q' is very close to the usual optimal Q, but is 
slightly larger since a' = a+ ee. In this example, 

Q' = 1.00129 . To interpret this result,recall that we 

are considering a uniform weighting of K at the various 
values of Q. K is relatively flat in the vicinity of its 
minimum, but increases somewhat more rapidly as Q decreases, 
than it does when Q increases, so we prefer to be "off 
optimal" on the high side. 

As a second example, suppose the contractor will 
usually produce and ship a quantity near the upper end of 
the allowed range. This is quite likely in real systems, 
since the manufacturer will generally "aim" for a quantity 


in the upper end of the range to maximize the sales profit. 
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Also, if the quantity actually produced is too small, he 
will make more, while if the quantity produced is too 
large, he supplies the maximum quantity allowed by the con- 
tract. This has a tendency to cause the delivery quantity 
distribution to be skewed to the right, with positive 
point probability at its right end point, even if the un- 
derlying production distribution is symmetric. For sim- 


plicity, assume a triangular distribution: 


£(Q) = a £(.99") = 0 , £(1.10') = yw. 
o* 2 


Now our problem is to minimize the expected value of K with 


respect to Q'. That is minimize 


1.1Q' 


K= JS (a/gtbg) 2007-90") ao 
-9Q' (Q') 
= 2, + b'Q! 
where a' = 10a - 45a 1n(1.222) = .98a, and b' = 10.035b - 


9.0b = 1.035b. This implies we should choose a smaller 
value for Q' than On Q' is approximately -9750.. Intui- 
tively, this is fairly easy to understand. We "expect" to 
get more than we order, so we "order short" to get the 
quantity we "really want". Note, however, that the changes 
above will be fairly nominal. The difference in the values 
of Q' and om is very small and the flatness of the cost 
function means that the cost saving is likely to be small 
so that it does not offset the extra cost of computation. 
The savings due to ordering Q' instead of Qu, in these two 


examples are .02 and .03 percent respectively. 
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We are now ready to consider the more difficult prob- 
lem of the cost of uncertainty. Suppose the lead time is 
long relative to the time between orders, so that several 
shipments are outstanding when we place an order. These é 
orders will be received during lead time. Normally, we 
would know the total asset picture exactly, but we now 
have an order quantity which is stochastic. In calculat- 
‘ing the optimal reorder point, we need to consider the 
difference between two random variables, instead of the 
difference between a constant and a random variable. The 
problem may be further complicated by the fact that one 
of these random variables is actually a sum of random 
variables which may not be independent. Assume for a - 
moment, however, that they are independent and define Y 
as the total quantity to be received during a lead time. 
The expected value of Y is n times the expected value of 
Q and it is relatively simple in theory to calculate the 
density of Y, as a convolution, and then from the density 
of X and this density of Y to calculate the density of 
Z=X - Y. We can use this density of Z in place of the 
density of demand X in the conventional models. Techniques 
for finding the densities of Y and Z can be found in 
Parzen [4]. 

However, it will not generally be possible to work 
with this sort of density, since the resulting function 
will be a special density and not a simple or tabulated 


one. Also it is not likely in practical systems that the 
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basic starting distributions or their parameters will be 
known. 

If we were fortunate enough to have adequate samvole 
data, we could estimate a distribution family and param- 
eters for the net inventory loss (or gain) during a lead 
time. Seldom will we have adequate data to do this, so it 
is usually necessary to assume a distribution, which seems 
appropriate, and estimate its parameters from the estimates 
of the parameters of f(x) and £(Q). 

It is reasonably clear that a fairly good estimate of 
the mean of net lead time demand will be the expected 
value of lead time demand minus the expected value of 
receipts (not n times the quantity ordered). Use of the 
parameter m = m 


x Q 
It is equally obvious that the variance * should be taken 


- nm, as the mean appears most reasonable. 


as greater than vie 


a problem, however. If we were convinced that the demand 


The question of how much greater is 


quantity during lead time and the delivery quantities QQ 


Qo sess Qn were each independent of the other, then 

o? = ‘* + no 9” would be appropriate. If we were convinced 
that the Q; were completely dependent, that is Q, =A for 
ald *®= lasses NM, then o* = Mee + ne On the other hand, 


the variance should be less than hy + no 2 if there is a 


Q 
tendency for things to even out. That is, if successive 
orders are negatively correlated, which may occur when the 


producer is trying to make corrections to his policy, the 


variance of total receipt quantity may be reduced below the 
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variance assuming independent delivery lot sizes. However, 
it appears reasonable to assume net leadtime demand has 
mean m =m, - nm. and variance 57 = 0 A + no 2 unless w 
ae Q ae Q % 
know the relationship between successive order sizes. 
The significant factor is that the variance of lead 
time demand should be increased by at least some amount and 


this will result in a higher safety level and/or greater 


risk of shortage, either of which increases operating costs. 
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CHAPTER 7 
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a OTHER RELAXATIONS OF THE CONVENTIONAL 
DELIVERY ASSUMPTIONS 
7.1 Stochastic and Variable Production Rates 
In the previous chapters it has been assumed that the 

production rate was either constant or a fixed function of 

time with no probabilistic element and that the inventory 

manager had no control over the rate of production. A 

production rate may be stochastic for several reasons. 

The production rate may be stochastic because it is a 

function of a stochastic process which generates imper- 

fect output. Machinery breakdown, strikes, emergency or- 

ders and other exogenous factors also influence the pro- 

~ duction rate. For example, the production rate (delivery 
rate) of a valve might be influenced by any of the follow- 
ing: the availability of brass for casting, timely delivery 
of brass to the foundry, other high priority foundry pro- 
jects, the actual percentage occurrence of flaws in the 
castings, variation in the quality of brass, foundry equipn- 
ment, or personnel performance, shipment from foundry to 
machine shop, machine breakdowns, high priority shov pro- 
jects, etc. 

t In this section we relax the assumptions in two ways. 
First, we will consider a stochastic production rate ¥ 
with probability density function f(¥), ¥ > 0. Consider 


a finite continuous production rate model for which both 
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demand and production are stochastic. See Figure 20. 
Let £(x,t) “be the probability density function of the 
demand X during an interval of length t. Let g(y,t) be the 
probability density function of the production Y during the 
interval t, derived from f£(¥). We assume all parameters of 
f(y,t) and q(y,t) are known. During the first part of a 
cycle, on hand stock builds up at the stochastic rate 
Z= Y - X. Once a total quantity of Q has been produced, 
production stops and the stock depletes at the rate X. 

Let h(z,t) be the density function of the net increase 


of inventory during the interval,t, defined as 


h(z,t) = f g(ztx,t) £(x,t) dx = Sf gly,t) f£ly-z,t) dy 


We use the product of the individual densities as the joint 
density since it is reasonable to assume X and Y are inde- 
pendent. It is quite realistic to allow the value of the 
random variable Z to be negative. This can easily happen 
during some time interval when demand is unusually large 
and production is unusually small. However, it is necessary 
for feasibility for Y to be stochastically greater than X, 
so that the expected value of Z2 = Y - X will be greater 
than zero. 

We can use the distribution h(z,t) to find the expected 
holding and shortage costs during the production phase of 
the cycle and the distribution f(x,y) to find the expected 


holding and shortage costs during the depletion phase of 


the cycle. However, we still have a problem, since both 





TY 


QUANTI 











FIGURE 20 


STOCHASTIC DEMAND AND PRODUCTION 
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the length of the production phase and the on hand quantity 
at the end of the production phase are stochastic. These 
values can not generally be found exactly, so we must use 
expected values for a3? Tar To, and on hand bei end of 


the production phase which is approximately r - AT, + 


0 
a . An exact model would be very complex and these 
approximations are probably quite adequate for some appli- 
cations. 

We can now set up an approximate model. It is exactly 
like the continuous production rate models of the literature, 
except that we use the distribution, h(z,t) during as and 
the usual distribution, f(x,t) during Ta only, instead of 
using f£(x,t) throughout the entire cycle with a constant 
ot daying T, “ 

Now consider a situation in which the production rate 
¥ can be controlled by the inventory manager. Specifically, 
we jecumes¥ = tec) where x is the cost per unit time which 
is passed on to the customer by the manufacturer. The cost 
of production as a function of the production rate is gen- 
erally accepted by economists to be convex from below, 
which implies some finite optimal rate greater than zero. 
There is usually an economy of scale effect at the lower 
production rates. Then at rates above the optimum, the 
production rate can be increased by the addition of ma- 
chinery, by subcontracting part of the work or by producing 
on overtime at premium wages. These costs can be expressed 


in terms of decreased or increased unit cost,but the effect 
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of this cost is clearer when isolated. Although a quick 
reference to the finite production rate model of section 
1.3 will suggest that a production rate of ¥ = A is opti- 

a mal, this applies only if we are overating on a portion 
of the ¥ range where the cost is near minimal. In general, 
to find an optimal we would add the cost C(¥) to the 


variable cost so that 


| K 


Ul 
QO 
: 
+ 

| 
A 


This can be solved for optimal values of ¥ and Q by the 
usual techniques. We also have a constraint to consider, 
since for feasibility, ¥Y* > \}. The same simple modifica- 
tion can be made to any finite production rate model. 
There is also a meaningful single period model to be 
é considered, when stock has dropped below the desired on 
hand level and we want to regain our steady state optimum. 
To do this, we wll have to increase ¥. For simplicity and 
clarity, we consider the model shown in Figure 21. We 
have b backorders and a deterministic demand rate 4. The 


“a 


backorder cost is 7 per unit per unit time. 





If we choose a ¥ that is much larger than dA we will 
save on the cost of backorders, but it will cost us more 
to achieve the rate ¥. If we let ¥ = A + € where € is 
. small, then it will cost us little extra to get the faster 
production, but the backorders will last for some time and 
be expensive. Is there an optimal ¥? 
The costs of importance in determining the optimal 


value of ¥ for this model are the costs of the backorders 
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FIGURE 21 


SINGLE PERIOD, VARIABLE ¥ MODEL 








during the time t until we are back in steady state opti- 
mal operation and the cost of the accelerated ¥ during the 
same period. Let cost = C(¥) be the cost of the rate ¥ 
(the inverse function of ¥ = g(x). 
Now the total cost to be minimized with respect to ¥ 
- a2 
K= C(¥)t + 5 bt = ate + ~— . 


For analytic convenience, we may make a change of variable 
to t= a and find the optimal value of t. If we are 


producing at or above the optimal rate ¥) ;then ,in this 
range,C(¥) is monotonically increasing in ¥. Typically, 
for some range up to ¥ max, we can increase our rate by 
some fixed amount for each hour of overtime. This results 


in a linear function. Suppose C(¥) = a+ a‘"¥. Then 


aon a a"h 

Cite) = Aa ah ot meer 
and 

4 Glee 4 Fy 

K = (ata')+ = Jt + 5 sit... 
Now 

aK , i = 

FFT ataaAr+aybdb=0, 


which is mot a function of t and implies t* = 0 or t* = ~, 


K = a'b + t(ata' A+5b) 


is monotonically increasing with t. Therefore t* = 0 and 


ye = dX + a= o, We interpret this to mean produce at ? 


max. 
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This is intuitively reasonable, since the linear cost 
of increasing production will imply a fixed cost of making 
up the shortage regardless of the time taken, and the faster 
the shortage is made up the lower the backorder cost. 

‘Now suppose C(¥) = ay’, This is a realistic approxi- 


mation of the cost curve above Yor the optimum production 


rate. Now 





‘ aby . le" 
~. WSR 2(¥-r) ’ 
and 
He» [2B - a ale: 
a¥ TAL ¢yaay 2 2(¥-a)2 
implies 
2% aE 


Similar results can be obtained for any function C(¥). 
7.2 Variable Lead Times and Expediting 

The advantages of shortened lead time are not apparent 
in deterministic models since there is no uncertainty in- 
volved. The major advantage of shortened lead time is due 
to the reduction in variance of lead time demand and is 
apparent in probabilistic models. However, in general, a 
reduced lead time has other benefits which are overlooked 


in most analytic models since the assumptions made exclude 


these effects. Reduced lead time réduces the risk of 





becoming overstocked and reduces the reaction time in 


situations where parameters are changing or subject to 
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change, and when stock is lost or destroyed. Reduced lead 
time has particular value when trends exist,or when sudden 
changes in demand rate or costs occur, or when the system 
is in trouble. 


fhe discussion of lead time control is divided into 


a 


* 
bese 


two parts. We will first consider routine reduction of 
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lead time as a continuing practice in a settied steady 
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state system with stochastic demand, This will be formu- 


lated as a trade off between the cost of decreasing lead 
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time and the expected cost saving'due to decreasing safety 
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level requirements. The second lead time reduction consid- 
ered will be referred to as expediting and differs from 
routine lead time reduction, since we will assume that an. 
improbable event has occurred during lead time and we face 
unusually large shortages prior to the planned delivery. 
This is formulated as a single period model. In real sys- 
tems the first model would be part of the routine stock re- 
plenishment computations. 

The second model is applicabie only when the system 
has gotten into temporary trouble due to unusually large 
demand, or loss of stock,or when making initial procure- 
ments on new items. To formulate a routine lead time op- 
timization model we can modify any of the convention mod- 
els in the literature by adding the decision variable t. 


We assume that a normal lead time, ty, is available at no 
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extra cost and that the cost of reducing lead time byt 

time units is c(t). We then replace the fixed value for 

lead time throughout the -model with T,.- t and add the cost 

_ term e(t). This results in a modified model to.be optimized : 
with respect to the original decision variables and t. 

For example, consider a very simple approximate mod- 
el as shown jin, Figure 22. The objective is to minimize the 
cost of. holding an inventory which is adequate to insure- 
that the probability of one or more shortages in a cycle is 
less than P, .If the normal lead time, The and the lead 
time demand distribution are known, we can find an optimal 
t as follows. The-holding cost is approximately |= + b | rc 
where b is the buffer stock required to keep the probability 5 
of shortage down to P. The cost of reducing lead time is 
c(t). Recall b will be. a function of lead time, lead time 
demand and P. The reorder-point, r, will be equal to b 
plus the expected lead time demand, At. The variable cost 
dependent upon t.is, 

aes s , K = c(t) + ICb(t). 

Suppose for example that c(t) = 200t7, the lead time 
demand is normal with mean At and variance ny AeSRZ0Y 
Ic = 40, P =..05 and Ty = o5. 

Now we know that to insure we have at most a .05° 
chance of shortage during lead time when the demand is nor- 
mally distributed,we must have stock on hand at the time. of 
order of at least thg mean lead time demand plus 1.65 stan- 


dard,deviations. .This means that b(t) is equal. to 1.65 
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FIGURE 22 


OPTIMAL LEAD TIME MODEL 





times the standard deviation Tay - t. We can now write 


K = 200t* + 40(1.65) (.5-t). 
Solving for the t which yields minimal costs, we find, 
£* = ,165. 

A single period model for expediting would be a fairly 
simple modification of the model above, but we would use 
the actual value of inventory position instead of the re- 
order point in the calculations. Allen and D'Esopo [5] 
present an analysis where an option exists to expedite 
delivery at a fixed extra cost. They determine an opti- 
mal policy for making a decision whether to expedite, but 
the analysis is limited by its assumptions. Generally, 
optimal inventory policies minimize long term exvected 
cost so there will often be specific points in time where 
improbable events have occurred and remedial action such 
as, expediting will result in additional savings. The fre- 
sence with which such situations arise in constrained bud- 
get systems, such as the military inventory systems,is quite 
sufficient to warrant further investigation. 
7.3 Multiple Source With Different Delivery Patterns 

Generally, an inventory manager has several sources 
for stock replenishment. It will be purely accidental if 
all sources offer the same delivery conditions. Some work 
has been done in considering multiple source procurement, 
notably the paper by W. J. Fabrycky and Jerry Banks [6] 
which considers a multiple source multiple item inventorv 


problem. Present practices usually choose the procurement 
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source solely on the basis of lowest item cost. It is 
ebyious from the considerations of the previous chapters 
that delivery conditions can influence total cost in a 
significant fashion. It may be cheaper over all‘ to.pav more 
for the items to get favorable delivery, particularly when 
the demand has a high variance to mean ratio or when the 
system is in short supply. Unfortunately, purchase rules 
generally preclude applying such enlightened management 
procedures. Dramatic changes. in bid evaluation | 
rules would be required. No attempt is made here to con- 
sider this facet of the inventory problem in detail. Sev- 
eral of the previous models can be applied in appropriate 
situations. At worst, it may be feasible to compute and 
compare the minimum cost per unit time for each available 
source of supply, and choose the source associated with 
the lowest of these minimums. 

Some comments on the Fabrycky-Banks paner are included 
here in order to provide an interface between the problem 
they approach and the problem of delivery variation. They 
consider the problem of multiple items with maltiple sources 
and give a dynamic programming technique for solution of the 
deterministic demand and lead time case. Their paper has an 
obvious application to the problem of relaxing delivery 
assumptions, although for each “source” the conventional 
delivery assumptions are made. A sepanebe order cost, unit 
price, leadtime and delivery rate are assigned to each 


“source”. “Source” is placed in quotation marks since 
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"source" in their formulation can as well apply to various 
modes of delivery.from a single manufacturer, and provides 


a method both for choosing between distinct sources and 


~ ‘between delivery options from a single source. It is di- 





rectly applicable to choosing optimally between specific 
alternative delivery schedules for which various lead times 
and delivery rates are available at different order costs 
and unit prices. 

Continuous models can be fit into their dynamic pro- 
gramming model by investigating a set of fixed options and 
then reinvestigating another set of more finely spaced op- 
tions in the vicinity of the minimum cost parameters. This 


procedure can be continued until the desired degree of 


accuracy is obtained. 
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CHAPTER 8 


SUMMARY AND CONCLUSIONS 

There is much work to be done in the area of relaxed 
delivery assumptions. Unfortunately, it does not seem 
likely that any general cohesive models will be forthcoming | 
in the near future. Unified treatment models become ex- 
ceedingly complex and intractable. The only feasible an- 
proach at present appears to be to generalize delivery as~ 
sumptions by setting up alternative models with different 
assumptions about the delivery conditions. 

Applicable models with more than one stochastic param- 
eter are always difficult to handle. Exact models can be 
formulated only when very restrictive assumeotions are made 
about the distributions of the random variables. Approxi- 
Mate solutions to models with simplifying assumptions made 
as required for algebraic and statistical tractability often . 
end up far afield. The major problem is that the ay ee 


are dictated by mathematical convenience and not by the na- 


“ture of the model. It may be much more difficult to assess 


the impact of such assumptions than to assess the impact of 
assumptions which are made about the model prior to any 
attempt at solution. 

Approximate models which are deterministic or incor~- 
porate only simple probabilistic considerations may pro- 
vide a manager with operation policies which are far supe- 


rior to the techniques previously used. 
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The basic ideas approached in this paper were that dé- 
liveries may be phased over time, that receipt (production) 
may he a stochastic process as well as issue (demand), 
that the anager may have control over certain parameters . 
usually taken as given (¥,t), and that the minimal cost 
policy may depend on the delivery conditions and options 
as well as the conditions usually Relea dtm eecount, 

There are several tentative conclusions that can be 
made based on the models and ideas in this paper. The 
conclusions drawn are: 

1. Phased delivery will usually result in larger or- 

der quantities and lower cost. 

2. The greater the inventory manager's control over 
the delivery conditions the cheaper a system can 
be run. 

Je Shortening lead time has a direct economic value 
in stochastic systems and is often worth paying 
for. 

4, When the manager has the option of determining the 
size of sellageeles as the contract progresses, the 
model will reduce to a standard delivery model 
with an adjusted cost of order. 

5.) Stechastve order quantities are unlikely to Sitect 
the optimal order policy significantly unless the 


distribution is very skewed or unless lead time is 


much larger than cycle time. 
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10. 


Unless constrained, the optimal contract periad 
will tend to be very large for most real a 
It will usually be a good approximation to place 

a contract for the longest period which common 
sense and the exogenous constraints allow. 

It is often possible to calculate the ontimal 
shipment size within a small error by using a 
modified purchase cost in the standard infinite 
delivery rate models. 

In terms of the defined parameters,optimal ship- 
ment lot size is an increasing function of B, the 
lot order cost, and is generally independent of A, 
the contract placement cost, while contract duca> 
tion is primarily a function of A and is relatively 
uninfluenced by B. 

When delivery conditioms vary among supplier's. 
bids, delivery conditions should he considered 

in choosing the “low-bid". The smallest unit 
price may not indicate the most economical source. 
Fractional delivery models,where fixed nercent- 
ages of the ordered quantity are delivered at 
fixed intervals, do not result in a change from 
the one lot delivery order quantity, if shortages 
are not allowed or are exceedingly -rare. 

In probabilistic phased delivery models each de- 
livery lot should he slightly Larger than the 


demand expected prior to the next delivery since 
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the increasing variance due to stochastic “demand 
(and delivery) will require deeper.safety stock 
to maintain comparable protection against stock 
out and achieve minimum cost. 


The subject for this paper was suggested by the 





author's experience in inventory management in the U. S. 
Navy Supply system. The Armed Services have become more 
cost conscious in recent years, and are among the leaders 
in computerized scientific inventory management. The scale 
of operations and the dollar volume makes even small per- 
-Gentage savings significant. There is a broad range of 
item characteristics represented and a wide selection of 
“delivery conditions. There is an obvious need for exten- 
sion of the conventional delivery assumptions to take ad- 
vantage of the unconventional delivery opportunities 
available. Of particular interest is the desirability of - 
annual contracts calling for phased delivery to keep the 
contracting cost down without raising on hand inventory 
levels. 
This paper represents a first assault on the problem. 
It is clear that the models and techniques proposed have 
considerable limitations. However, they provide a start- 
ing point from which to proceed. It is also hoped that 


some provide a means to determine whether the anticipated 


savings will justify the expense of detailed analysis. 
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